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vs.  Dimensionless  Film  Thickness  (Q  =  0. 175W, 

So  =  3.6  nm,  Kq  =  2.84  x  102  N/m2), 

Figure  5.17  Plot  of  Dimensionless  Pressure  and  Curvature  91 
vs.  Dimensionless  Film  Thickness  (Q  *  0.54W, 

So  *  2.9  nm,  tto  =  5.35  x  102  N/m2). 

Figure  5.18  Plot  of  Dimensionless  Pressure  and  Curvature  91 

vs.  Dimensionless  Film  Thickness  (Q  =  0.77W, 

So  =  2.7  nm,  Kq  =  1.85  x  10^  N/r 

Figure  5.19  Plot  of  Dimensionless  Pressure  id  Curvature  9  2 

vs.  Dimensionless  Film  Thickness  (Q  =  O.OW, 

So  =  17.2  nm,  JCo  =  3.88  x  lOO  N/m2). 

Figure  5.20  Plot  of  Dimensionless  Pressure  and  Curvature  9  2 
vs.  Dimensionless  Film  Thickness  (Q  =  0.05W, 

So  *  5.9  nm,  «o  =  7.53  x  10*  N/m2). 

Figure  5.21  Plot  of  Dimensionless  Pressure  and  Curvature  9  3 

vs.  Dimensionless  Film  Thickness  (Q  *  0.075W, 

So  =  5.4  nm,  ?to  =  1  28  x  102  N/m2). 

Figure  5.22  Plot  of  Dimensionless  Pressure  and  Curvature  9  3 

vs.  Dimensionless  Film  Thickness  (Q  »  O.IW, 

5q  3  4.9  nm,  Jto  •  2.42  x  102  N/m2). 
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Figure  5.23  Plot  of  Dimensionless  Pressure  and  Curvature  9  4 
vs.  Dimensionless  Film  Thickness  (Q  ■  0.125W. 

6o  -  4.5  nm,  Kq  «  3.68  x  102  N/m2). 

Figure  5.24  96 

Dimensionless  Mass  Flux,  vs.  Dimensionless  Film  Thickness.  ti. 

Case  1  Q  *  O.OW,  8o  »  15  nm.  m'd  «  6.83  x  10-5  kg/m2»sec' 

^Case  2  Q  »  O.IW.  5o  -  6.4  nm.  m'4  «  8.50  x  10-4  kg/m2.sec^ 

Figure  5.25  96 

Dimensionless  Mass  Flux.  ^  vs.  Dimensionless  Film  Thickness,  n. 

Case  1  Q  »  O.IW,  6o  *  6.4  nm.  m‘4  a  8.50  x  10*4  kg/m2»sec 

Case  2  Q  »  0.3W.  6o  =  5.9  nm,  m‘4  »  9.99  x  10*4  kg/m2»sec^ 

Figure  5.26  9  7 

Dimensionless  Mass  Flux,  vs.  Dimensionless  Film  Thickness.  t|. 

'case  1  Q  -  0.2W.  6o  »  6.2  nm.  m‘d  «  6.63  x  10*4  kg/m2.sec  ' 
^Case  2  Q  *  0.5W,  «o  *  5.65  nm.  m‘4  »  1.15  x  10*3  kg/m2.sec^ 

Figure  5.27  9  7 

Dimensionless  Mass  Flux.  Nf  vs.  Dimensionless  Film  Thickness,  r|. 

Case  1  Q  »  O.OW,  6©  =  11.5  nm,  m'*!  =  1.92  x  10*3  kg/m2»sec 
Case  2  Q  »  0.35W,  Sg  »  4.5  nm,  m***  =  1.04  x  10*2  kg/m2»secj 

Figure  5.28  9  8 

Dimensionless  Mass  Flux,  M  vs.  Dimensionless  Film  Thickness,  t\. 

Case  I  Q  ■  0.25W,  5©  »  4.8  nm,  m*4  »  i.27  x  10*2  kg/m2.sec 
^Case  2  Q  =  0.5W,  fig  *  4.0  nm,  m'^l  »  2.81  x  10*2  kg/m2»sec  ^ 

Figure  5.29  98 

Dimensionless  Mass  Flux,  M  vs.  Dimensionless  Film  Thickness,  n- 
Case  1  Q  «  O.OW,  8g  ■  13.6  nm,  m'4  »  5.64  x  10*®  kg/m2»sec 
^Case  2  Q  •  O.IW,  8g  ■  4.7  nm,  m'®  ■  4.73  x  10*5  kg/m2»sec  ^ 
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Figure  5.30  9  9 

Dimensionless  Mass  Flux,  vs.  Dimensionless  Film  Thickness,  n 

Case  1  Q  =  0.075W,  So  »  5.2  nm,  m'*^  =  3.84  x  lO'^  kg/m^'sec^ 

Case  2  Q  »  0.125W,  So  =  4.4  nm,  m‘4  s  5.39  x  10-5  kg/m^.sec^ 

Figure  5.31  99 

Dimensionless  Mass  Flux,  ^  vs.  Dimensionless  Film  Thickness,  ii. 

Case  1  Q  =  O.OW,  So  =  17.2  nm,  m‘4  s  6.24  x  10*"^  kg/m2*sec 
^  Case  2  Q  =  O.IW,  So  =  4.9  nm,  m*4  =  4.18  x  10‘5  kg/m7*sec  ^ 

Figure  5.32  100 

Dimensionless  Mass  Flux.  ^  vs.  Dimensionless  Film  Thickness,  ii. 

Case  1  Q  *  0.05W.  So  =  5.9  nm,  m‘4  =  1.30  x  10-5  kg/m2»sec' 
Case  2  Q  *  O.IW,  So  =  4.9  nm,  m*4  a  4.18  x  10*5  kg/m2*sec  ^ 

Figure  5.33  100 

Dimensionless  Mass  Flux,  ^  vs.  Dimensionless  Film  Thickness,  t]. 

Case  I  Q  =  0.075W,  So  *  5.4  nm,  m>4  »  2.21  x  10*5  kg/m2«sec' 

^Case  2  Q  =  0.125W,  So  =  4.5  nm.  m‘4  =  6.64  x  10*5  kg/m2*sec^ 


Figure  5.34 

Integral  Evaporative  Heat  Sink  vs.  Film  Thickness 
(Q  »  0.5W  for  all  the  cases). 
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SECTION  1 


SUMMARY 

The  long-range  objectives  of  this  research  program  are  to 
identify,  evaluate,  and  correlate  important  characteristics  of  the  heat 
transfer  process  in  evaporating  thin  liquid  films  that  can  be 
subsequently  used  to  optimize  their  performance  as  heat  sinks. 
Herein,  a  thin  liquid  film  in  the  shape  of  a  meniscus  is  emphasized. 

In  the  macroscopic  experimental  pan  of  the  study,  a  relatively  large 
extended  meniscus  was  studied  at  low  evaporation  rates.  However, 
we  believe  that  the  characteristics  obtained  for  low  evaporation 
rates  and  relatively  large  systems  should  also  represent  the 
phenomena  in  smaller  systems  with  higher  evaporation  rates.  In 
the  microscopic  experimental  and  theoretical  pan  of  the  study,  the 
contact  line  region  of  the  meniscus,  which  is  a  small  (in  general, 
thickness  less  than  the  order  of  S  p)  but  critical  part  of  the  extended 
meniscus,  was  studied.  As  the  meniscus  decreases  in  size,  the 
contact  line  region  would  represent  a  larger  ponion  of  the  total  heat 
sink.  In  the  latter  case,  the  stability  of  the  system  would  increase 
because  of  the  closer  proximity  of  the  solid  surface  which  (along 
with  the  liquid  properties)  passively  controls  the  shape  of  the 
meniscus.  Therefore,  the  experimental  approach  is  to  study  the 
characteristics  of  the  evaporating  meniscus  in  a  sufficiently  large 
system  that  allows  the  measurements  and  subsequent  modeling  of 
the  results.  Using  the  kinetic  theory  ideal  interfacial  heat  flux  as  a 
reference,  the  relative  effects  of  temperature,  pressure,  and 
conduction  on  the  interfacial  heat  flux  are  determined.  Having 
obtained  an  experimental  confirmation  of  the  modeling,  the  model 
can  then  be  used  to  predict  the  performance  of  smaller  size  but 
higher  flux  systems.  We  find  that  the  characteristics  of  an 
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evaporating  thin  film  which  are  controlled  by  interfacial  phenomena 
are  extremely  complicated  but  tractable. 

Initially,  we  proposed  to  emphasize  the  "interline  heat  flow 
parameter."  Because  of  the  theoretical  results  of  the  research,  this 
objective  has  been  expanded  to  include  the  "capillary  pumping 
parameter"  for  completeness.  The  interline  heat  flow  parameter  is 
important  in  ultra-thin  films  where  the  liquid-solid  interfacial  force 
field  predominates,  whereas  the  "capillary  pumping  parameter"  is 
important  in  thicker  films.  In  the  transition  region,  they  are  both 
important.  In  the  theoretical  modeling  of  the  complete  meniscus,  we 
then  found  it  convenient  to  replace  these  two  dimensional 
parameters  with  a  single  dimensionless  group,  e.  Although  the 
analysis  is  restricted  by  the  lubrication  approximation  to  the  small 
slope  region  of  the  meniscus,  this  dimensionless  group  has  many 
facets  which  enhance  the  characterization  of  the  entire  meniscus. 

The  liquid-solid  interfacial  force  field  and  temperature  difference  fix 
the  characteristic  contact  line  thickness,  8^,  which  is  of  critical 
imponance  because  it  controls  the  rest  of  the  meniscus.  The 
measured  profile  gives  8^  and  e.  Finally,  we  added  the  effect  of 
thermal  conduction  in  the  liquid  to  the  model. 


The  integral  evaporative  heat  sink  for  the  region  where  the 
thickness  is  less  than  8.  Qg,  can  be  obtained  by  evaluating  Eq.  (1.1)  at  8. 


Q8  = 


8  3K’ 

f-  ^ 

AAh^ 

1  J 

3 

1  J 

§1 

8 


where 


(l.l) 


-  is  the  capillary  pumping  parameter. 


AAhm 

-  is  the  interline  heat  flow  parameter, 
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The  corresponding  dimensionless  equation  for  the  integral  heat  sink 
evaluated  at  i\  is 


Q 

where 


*  3  d^ 


(1.2) 


Q=  -Qsxov/Mhm,  ♦  =  -(n*^  +  e^").  e  =  (o  5q  V  b)/(-A), 

ti=5/5o.  §  =  x// 

Equation  (1.1)  demonstrates  that,  if  the  physical  properties  of  the 
system  are  known,  the  experimentally  measured  film  thickness 
profile  can  be  used  to  calculate  the  integral  evaporative  heat  sink. 
Therefore,  one  of  the  major  experimental  objectives  of  the  research 
was  to  determine  the  film  thickness  profile  as  a  function  of  the 
experimentally  varied  integral  evaporation  rate.  In  addition  the 

unknown  modified  Hamaker  constant.  A,  had  to  be  determined  in 
situ  because  it  is  a  strong  function  of  the  surface  history  of  the 
system.  We  note  that,  although  the  Hamaker  constant  is  a  basic 
constant  which  can  be  theoretically  calculated  for  ideal  surfaces, 
most  engineering  surfaces  are  not  ideal.  Therefore  the  Hamaker 
constant  has  to  be  experimentally  obtained. 

There  are  two  main  aspects  of  the  research:  experimental  and 
theoretical.  We  will  present  the  experimental  results  in  Section  3. 
This  will  be  followed  by  a  presentation  of  the  theoretical  results  in 
Section  4  where  the  data  are  successfully  used  to  evaluate  the 
theoretical  model.  A  primer  on  surface  physics  which  is  germane  to 
the  research  and  an  example  demonstrating  the  use  of  the  results  of 
the  modeling  is  given  in  the  Appendices  A  and  B.  Because  of 
experimental  necessities,  there  are  also  microscopic  and  macroscopic 
aspects  of  the  research. 
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1.1  Experimental 


At  the  microscopic  level,  using  ellipsometry  and 
microcomputer  enhanced  video  microscopy  (based  on 
interferometry),  the  film  thickness  profile  was  measured  as  a 
function  of  the  evaporation  rate  for  polar  and  nonpolar  fluids.  These 
measurements  were  made  in  the  contact  line  region  where  the 
thickness  was  less  than  3  p.  Therefore,  the  integral  evaporative  heat 
sink,  Qg,  could  be  obtained  for  this  region  using  Eq.  (1.1).  We  note 

that  the  imponant  characteristic  thickness  at  the  contact  line,  5^, 

was  experimentally  measured.  We  also  note  that  a  portion  of  the 
line  tension  at  the  contact  line  (force/length),  which  keeps  the 
meniscus  from  dewetting,  is  given  by  ^5^^.  Experimental 

limitations  did  not  allow  the  full  potential  of  the  evaporating 
meniscus  to  be  measured  since  the  smallest  value  of  5^  achieved 

was  2.7  nm,  which  is  substantially  larger  than  the  presumed  limit  of 
a  monolayer. 

At  the  macroscopic  level,  a  heat  balance  was  obtained  by 
comparing  the  total  integral  evaporation  rate  for  the  entire  meniscus 
based  on  the  integral  input  mass  flow  rate,  Qj,  with  that  based  on 

the  temperature  profile.  The  above  was  achieved  using  a  circular 
heat  transfer  cell  with  a  capillary  feeder  which  was  designed  and 
developed  for  this  project. 

1.2  Theoretical 


As  described  in  Section  4,  the  following  theoretical  equation 
was  derived  to  model  the  interfacial  mass  flux  associated  with  fluid 
flow  in  an  evaporating  thin  film. 


I  A  f  3 
3  d5  I’’  d^J 


1 

1  +  KTl 


(1  +  4>) 


4 


(1.3) 


where  the  parameter  k  is 


aAhm  Sq 
^  *  k 


This  fourth  order  nonlinear  differential  equation  for  the 
dimensionless  thickness  has  two  parameters,  c  and  k.  The 
parameter  e  represents  the  effects  of  intermolecular  forces  and  k 
represents  the  resistance  to  thermal  conduction  in  the  meniscus. 

The  only  unknown  in  the  parameter  k  was  obtained  by  measuring 
5q.  The  experimentally  measured  thin  film  thickness  profile  data 

was  successfully  fitted  with  the  model  equation  to  obtain  the 
important  dimensionless  parameter,  e,  which  includes  both  the 
surface  tension  and  the  Hamaker  constant.  This  number  combines 
the  properties  of  the  capillary  pumping  parameter  and  the  interline 
heat  flow  parameter.  Knowing  the  classical  properties  of  the  system, 
the  Hamaker  constant  is  obtained  from  e.  Once  the  parameter  e  and 
the  contact  line  film  thickness,  5^,  are  known,  the  ideal  heat  flux, 

q*^,  and  the  liquid-vapor  and  solid-vapor  interfacial  temperature 
differences  can  be  calculated.  The  interfacial  temperature 
difference  cannot  be  measured  directly  because  of  the  exceptionally 
small  thickness  of  the  liquid-vapor  interface  across  which  the 
density  changes  from  that  of  the  liquid  to  that  of  the  vapor  and 
because  the  temperature  difference  is  extremely  small. 


The  right-hand  side  of  Eq.(1.3)  is  both  the  dimensionless 
interfacial  mass  flux,  I^,  and  the  dimensionless  interfacial  heat  flux: 


M 


m 


1 


«id  1  +  xq 


(1  +4>) 


m 


(1.4) 


The  ideal  interfacial  mass  flux  based  on  kinetic  theory,  m‘4,  is  a 
convenient  reference  flux.  We  find  that,  because  the  interfacial 
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force  field  (represented  by  ^  which  has  a  negative  value)  lowers  the 
vapor  pressure  and  conduction  in  the  liquid  (represented  by  Kq) 
lowers  the  temperature  at  the  liquid-vapor  interface,  the  ideal  heat 
flux  based  on  kinetic  theory  and  the  solid-vapor  interfacial 
temperature  difference,  AT,  is  greater  than  the  maximum  heat  flux 
in  the  evaporating  meniscus.  We  note  that  the  results  demonstrate 
that  these  effects  are  large.  The  results  also  indicate  that  high 
interfacial  mass  (heat)  fluxes  can  be  obtained  in  small  systems 
(small  KT>)  if  they  are  stable.  We  note  that  stability  should  increase 
in  small  wetting  systems  because  the  relative  importance  of  the 
solid  substrate  increases.  However,  this  complex  question  has  not 
been  addressed  and  the  experimental  system  was  relatively  large. 

The  resulting  theoretical  flt  of  the  data  is  also  used  to 
obtain  the  experimental  integral  heat  sink,  Qg,  and  the  heat  flux 

distribution.  However,  since  the  model  differential  equation  was 
derived  using  the  lubrication  approximation,  it  can  only  be 
accurately  used  in  the  region  where  (5’)^  «  1.  This  region  represents 
only  a  portion  of  the  complete  macroscopic  evaporating  meniscus. 

As  a  result  of  the  difference  in  areas,  a  gap  remains  between  the 
macroscopically  measured  Qj  and  the  microscopically  obtained  Qg  . 

We  note  that  this  does  not  detract  from  the  main  conclusion  of  the 
research:  the  experimental  results  are  in  agreement  with  the 
proposed  model  of  the  evaporating  meniscus.  Although  considerable 
additional  research  needs  to  be  accomplished  concerning  various 
fluids,  surface  preparation,  smaller  constrained  systems,  stability, 
and  higher  heat  fluxes,  the  basic  characteristics  of  the  evaporating 
meniscus  have  been  experimentally  measured  and  theoretically 
modeled. 
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1.3  Application  and  Surface  Physics 

In  the  Appendix  B,  an  example  of  the  use  of  some  of  the 
model  equations  to  evaluate  the  effect  of  the  contact  line  heat  sink 
at  the  exit  of  a  capillary  tube  on  the  capillary  pressure  drop  for  fluid 
flow  is  presented.  Only  the  effect  of  the  disjoining  portion  of  the 
evaporating  meniscus  on  the  capillary  pressure  is  evaluated  in  the 
application  section.  We  note  that  the  dimensionless  group  e  is 
needed  to  relate  the  size  of  the  heat  sink  in  the  example  with  that  in 
the  experimental  ponion  of  the  research.  The  importance  of  the 
parameters  e  and  k  is  demonstrated. 

Since  the  use  of  the  Hamaker  constant  to  model  evaporating 
ultra-thin  films  is  relatively  new,  a  primer  on  surface  physics  which 
is  germane  to  the  research  is  also  included  in  the  Appendix  A. 

1.4  Major  Conclusions 

I)  The  experimental  results  obtained  for  the  relatively  low 
heat  fluxes  studied  generally  confirm  the  basic  model  presented. 
Therefore,  the  model  can  be  used  to  predict  the  performance  of 
other  systems  and/or  operating  conditions. 

II)  We  find  that  there  are  significant  resistances  to  heat 
transfer  in  a  small  meniscus  due  to  interfacial  forces,  viscous 
stresses,  and  thermal  conduction.  Therefore,  the  ideal  heat  flux 
based  on  kinetics  theory  cannot  be  obtained. 
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SECTION  2 


INTRODUCTION 

When  a  completely  wetting  liquid  film  spreads  on  an 
inclined  solid  substrate  a  "contact  line"  region  is  formed  where  the 
vapor,  liquid  and  solid  (with  an  adsorbed  film)  phases  are  in  close 
proximity.  At  equilibrium  on  a  horizontal  surface  a  uniform  film  of 
constant  thickness  would  form.  Alternatively,  a  partially  spreading 
fluid  would  form  an  apparent  finite  contact  angle.  We  currently 
focus  on  completely  wetting  systems  for  experimental  and 
theoretical  advantages,  and  their  basic  importance.  The  transport 
processes  and  interfacial  phenomena  occurring  in  the  microscopic 
contact  line  region  (thickness  5  <  10*5  m]  are  controlled  by  surface 
forces  and  involve  a  delicate  balance  of  temperature  and  of 
capillary,  dispersion  and  hydrostatic  forces.  On  the  other  hand, 
large  stresses  are  present  in  this  interfacial  region  which  includes 
the  varying  force  fields  associated  with  both  the  liquid-vapor  and 
the  liquid-solid  interfaces.  The  stress  field  in  the  contact  line  region 
controls  the  rest  of  the  meniscus.  The  technical  applications  of  this 
generic  region  are  legion. 

2.1  Literature  review 

2.1.1  General 

Although  considerable  theoretical  progress  [e.g.,1-26]  and 
exciting  experimental  developments  [e.g., 27-42]  have  been  reported, 
a  complete  and  consistent  understanding  of  the  complex 
physicochemical  phenomena  in  the  contact  line  region  is  still  to  be 
achieved.  The  primary  obstacles  to  this  endeavor  are  the  smallness 
of  the  region  under  study,  the  large  stress  gradients,  and  its  extreme 
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sensitivity  to  the  thickness  profile,  the  temperature  gradient  and  to 
contamination.  Although  these  factors  are  difficult  to  include  in 
theoretical  analyses  and  vary  greatly  from  one  experimental  study 
to  another,  a  general  understanding  of  the  contact  line  is  evolving. 

Consider  Fig.  2.1  which  represents  an  inclined  flat  plate 
immersed  in  a  liquid  at  an  angle  6  to  the  horizontal.  The  plate  is  in 
a  closed  cell  and  a  spreading  liquid  in  equilibrium  with  its  own 
vapor  forms  a  zero  contact  angle  with  the  plate.  Except  for  the 
inclined  support  plate,  the  liquid  is  not  constrained  by  solid  surfaces 
and  the  liquid-vapor  interface  is  free  to  adjust  to  the  "pressure 
jump"  at  the  liquid-vapor  interface.  Therefore,  the  thickness  profile, 
5,  is  a  function  of  the  pre«sure  jump  at  the  interface  and  a 
measurement  of  this  profile  yives  the  pressure  field.  At 
equilibrium,  the  pressure  field  is  hydrostatic.  With  evaporation,  the 
change  in  the  profile  causes  a  change  in  the  pressure  field  that,  in 
turn,  causes  fluid  flow  to  replenish  that  evaporated.  Since  the 
pressure  field  depends  on  the  intermolecular  force  field,  extremely 
large  pressure  gradients  are  possible  without  cavitation.  A  steady 
state  evaporating  meniscus  has  been  experimentally  observed  [e.g., 
32,  33].  At  high  evaporative  fluxes,  an  oscillating  meniscus  has  also 
been  observed.  At  very  high  fluxes,  the  meniscus  becomes  unstable. 

We  note  that,  due  to  the  concept  of  internal  pressure, 
the  concept  of  pressure  jump  at  the  liquid-vapor  interface  is  not 
straightforward.  In  addition  the  internal  pressure  is  modified  by 
the  presence  of  the  solid  substrate.  However,  as  discussed  below, 
models  have  been  developed  and  experimentally  verified  for  the 
pressure  jump  under  some  conditions  [e.g. ,1,2,8,9,18,22, 27, 32,36]. 
Past  theoretical  [e.g. ,1,5, 10,1 1,16,36]  and  experimental 
[e.g. ,36,40,4 1,42]  work  demonstrated  that  there  is  a  narrow 
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Figure  2.1  An  Inclined  Flat  Plate  Immersed  in  a  Liquid 
at  an  Angle  9  to  the  Horizontal. 
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transition  region  between  the  region  dominated  by  capillary 
pressure  and  that  dominated  by  disjoining  pressure. 


2.1.2  Equilibrium  Systems  (Theoretical) 

Mechanical  equilibrium  requires  a  balance  of  dispersion, 
capillary  and  hydrostatic  forces.  Hence  the  following  form  of  the 
chemical  potential  (per  unit  volume)  is  zero  at  the  liquid-vapor 
interface  [e.g.,  8,19,22,36]. 


n  +  oK  -  pjgH  =  0 

(2.1) 

AP  (Pressure  Jump  at  Liquid-Vapor  Interface)  = 

Pi  -  Pv  =  -  oK  -  n 

(2.2) 

d25  /'dSV  *3/2 

K  (Curvature)  =  ^  [1  +  1 

(2.3) 

n  (Disjoining  Pressure)  = 

(2.4) 

In  these  equations,  5  represents  the  film  thickness,  p/  represents  the 
density,  a  represents  the  surface  tension  and  B  a  modified  Hamaker 
constant  (  B<0  for  completely  wetting  systems).  In  the  limit  of  a 
very  thin  film  of  a  pure  simple  fluid  (nonretarded  region)  n=3  and 
B=A/6tc  =  A,  in  which  A  is  the  classical  Hamaker  constant;  while  in 
the  thick  film  region  (retarded  region)  n=4  and  B  is  a  dispersion 
constant.  Calculations  of  the  van  der  Waals  interaction  between 
different  materials  show  that  these  are  approximate  models  since 
the  constants  are  weak  functions  of  the  film  thickness  [e.g.,2,36]. 

The  effect  of  buoyancy  is  neglected  as  experimental  conditions  are 
far  from  the  critical  point  (p|  »  Py).  Knowing  the  surface  tension,  a 

measurement  of  the  thickness  profile  gives  the  constant  B  and, 
therefore,  the  characteristics  of  the  system  at  the  start  of  an 


experiment  [36].  For  simple  clean  ideal  systems  the  measured  value 
of  the  Hamaker  constant  can  be  compared  with  the  theoretical  value 
calculated  using  the  frequency  dependent  dielectric  susceptibilities 
of  the  three  phases.  For  more  complicated  real  systems,  the  value  of 
the  Hamaker  constant  can  be  determined  in  a  closed  cell  in  situ  at 
the  start  of  an  experiment. 

We  note  that,  if  the  extended  meniscus  is  formed  at  the  exit  of 
a  small  capillary,  the  pressure  jump  can  be  extremely  large.  In  the 
capillary  case,  the  constraint  of  the  solid  can  be  used  to  reduce  the 
size  of  the  system  and  increase  the  relative  importance  of  the 
interfacial  force  terms.  This  would  increase  the  stability  of  the 
meniscus  and  the  attainable  heat  flux.  Because  of  the  current 
experimental  state  of  the  art,  a  relatively  large  meniscus  at  the  exit 
of  a  circular  gap  was  studied  herein  using  a  newly  designed  circular 
capillary  feeder. 

2.1.3  Nonequilibrium  Systems  (Theoretical) 

For  general  situations.  Equations  (2.6-2.10)  have  been  derived 
using  the  lubrication  approximation  to  the  one-dimensional  Navier- 
Stokes  equations  while  ignoring  surface  tension  gradients  [e.g., 
8,9,11,18,19,24,26,33].  A  more  complete  description  and  derivation 
of  the  equations  used  herein  are  presented  in  Section  4. 

r  (Mass  Flow  Rate  per  Unit  Width)  = 

5  -53  d 

Pi  37  57  IP’’*' 

0 

APg  =  AP  +  p,gx  sine  +  p,g(8-d)  cos0  (2.6) 

d  =  Film  thickness  at  pool,  where  x  =  0. 


In  the  following  equation  for  the  conservation  of  mass,  m 
represents  the  interfacial  mass  flux  and  t  represents  the  time. 


(2.7) 

(2.8) 


Using  kinetic  theory  and  interfacial  thermodynamics  with 
combined  Clapeyron  and  Kelvin  equations,  the  mass  flux  of  vapor 
leaving  the  liquid  vapor  interface,  m,  was  modeled  as  a  function  of 
the  temperature  and  pressure  jumps  at  the  interface  in  [e.g.,8,9, 
11,19,44].  The  resulting  equation  can  be  called  either  an  extended 
Kelvin  equation  or  an  extended  Clapeyron  equation. 

m  =  a(  -  Ty  )  +  b(APg)  (2.9) 


b 


Cl 


f  ^  1 

1/2 

/Ah^ 

(^2jiRTiv^ 

[rTiv  j  ’  a  -  I 

(2.10) 


where  Ci  is  the  evaporation  coefficient.  Ah  is  the  enthalpy  of 
vaporization  per  unit  volume,  M  is  the  molecular  weight,  is  the 

reference  vapor  pressure,  V|  is  the  molar  volume  of  the  liquid  and 
(T/v  -  Tv)  is  the  characteristic  temperature  jump  at  the  liquid- vapor 
interface.  These  equations  demonstrate  that  the  processes  of  fluid 
flow  and  phase  change  both  depend  on  the  intermolecular  force 
field  and  are  strongly  coupled.  We  note  that  the  extended  effusion 
equation,  Eq.  (2.9),  predicts  that  extremely  large  evaporative  heat 
fluxes  are  possible  if  the  system  does  not  cavitate.  Very  little  is 
known  about  this  generic  area.  We  also  note  that,  as  the  system 
shrinks  in  size  with  evaporation,  extremely  large  shear  stresses 
occur  because  of  flow  in  very  thin  films.  Therefore,  in  order  to 


obtain  a  complete  understanding  of  the  transport  processes  in  an 
evaporating  extended  meniscus,  it  is  important  to  evaluate 
equilibrium  systems  (to  calibrate  the  system  at  the  start  of  the 
experiments  and  evaluate  the  interfacial  force  field),  near¬ 
equilibrium  systems  (to  evaluate  models  at  low  fluxes  where  the 
system  does  not  oscillate)  and  high  flux  systems  (to  evaluate  the 
system  when  the  system  oscillates  and  cavitation  occurs).  The  first 
two  are  addressed  herein. 

An  examination  of  the  relevant  physical  properties  of  many 
fluids  shows  that  these  systems  are  extremely  sensitive  to  a  small 
AT.  On  the  other  hand,  we  can  theoretically  demonstrate  that  a 
small  change  in  AT  is  equivalent  to  a  large  change  in  the  stress  field. 
Herein  we  are  primarily  concerned  with  spreading  systems.  Thus,  fl 
and  oK  are  both  positive  and  the  vapor  pressure  of  these  films  is 
less  than  that  of  a  bulk  liquid  at  the  same  temperature.  The  above 
equations  model  the  evaporative  flux  as  a  function  of  both 
temperature  and  pressure  (film  thickness  and  curvature).  Derjaguin 
and  Zorin  [27]  experimentally  demonstrated  that  thin  superheated 
(relative  to  a  bulk  liquid  phase)  flat  adsorbed  liquid  films  of  a 
particular  thickness  on  a  solid  substrate  can  exist  in  equilibrium  (fn 
=  0)  with  the  saturated  vapor  of  the  bulk  liquid  phase.  Thicker 
films  at  the  elevated  temperature  would  have  a  higher  vapor 
pressure  and  therefore  evaporate  until  they  reach  the  equilibrium 
thickness  at  which  the  additional  liquid-solid  interfacial  force  field 
reduces  the  vapor  pressure.  Using  Eq.  (2.9)  for  an  adsorbed  film  on 
a  flat  solid  surface  above  in  the  vapor  space  but  not  in  contact  with 
the  bulk  liquid  (  K  =  0  and  H  ^0)  we  find  this  thickness  to  be 


We  note  that,  in  this  region,  AT  *  T/v  -  Ty  =  T*  -  Ty.  This  thickness  is 

also  used  to  define  the  characteristic  pressure  difference  (disjoining 
pressure  at  n^). 

„  -B  Ah  AT 
~  T 

*0  (2.12) 

Using  an  example  from  Ref.  [19],  the  equivalent  effect  of  a  AT  =  0.1  K 
for  Octane  on  Si02  at  298  K  with  B  s  A/6jt  =  -4.35  x  10'22  J,  n  =  3,  is 
given  by  Flo  =  8.5  x  lO^  N/m2  (5o  =  1.72  nm,  H  =  0);  K©  =  4  x  10®  m 
(large  5^,  H  =  0);  or  H  =  12.4  m  (  K=0,AT  =  0).  It  is  evident  from  these 

preliminary  calculations  and  confirmed  herein  that  small 
temperature  differences  cause  extremely  thin  films  and/or  large 
curvatures  and  large  pressure  gradients.  The  thickness  is 

important  because  this  is  the  boundary  condition  for  an  evaporating 
meniscus  at  the  contact  line  [e.g.,  9,  19].  The  upper  limit  on  AT  for 
this  boundary  condition  is  unknown.  However,  it  can  be  sufficiently 
large  to  allow  very  high  heat  fluxes.  We  note  that  Ref.[19]  was 
written  under  the  auspices  of  this  contract.  Herein,  we  use 
experimental  data  to  confirm  the  characteristics  predicted  by  this 
model. 

2.1.4  Equilibrium  Systems  (Experimental) 

Essentially,  three  different  experimental  techniques  have  been 
used  to  study  the  equilibrium  disjoining  pressure  of  a  simple  system 
like  an  alkane  on  a  wetting  substrate:  I.)  Derjaguin  and  Zorin  [27], 
Zorin  et  al.  [34]  and  Gee  et  al.  [37]  studied  the  adsorption  of  alkanes 
on  substrates  at  high  disjoining  pressures  (i.e.,  very  thin  adsorbed 
films  which  were  thinner  than  those  studied  herein).  Derjaguin  and 
Zorin  varied  the  substrate  temperature  whereas  Gee  et  al.  varied  the 
vapor  pressure;  II.)  Blake  [30]  measured  the  disjoining  pressure  of 


alkanes  on  a-alumina  by  pressing  a  gas  bubble  against  the  substrate 
at  low  disjoining  pressures  which  are  not  easily  attainable  by  vapor 
adsorption;  III.)  Ingram  [29]  studied  very  low  disjoining  pressures 
by  measuring  the  thicknesses  of  alkane  films  which  remained  on  a 
silica  plate  after  draining.  Truong  and  Wayner  [36]  and  Sujanani 
and  Wayner  [40]  used  the  set-up  presented  in  Figure  2.1  to 
measure  the  profile  in  both  the  low  disjoining  pressure  region  just 
above  the  transition  region  and  the  capillary  pressure  region.  The 
results  presented  in  [27,34,36,37]  agree  that  the  DLP  theory  [2] 
correctly  predicts  the  value  of  the  Hamaker  constant.  Ingram  found 
that  octane  and  higher  alkanes  gave  much  thinner  films  after 
drainage  than  predicted  whereas  the  lower  alkanes  gave  the  thicker 
films  expected  [29].  He  attributed  this  discrepancy  to  either 
contamination  (possibly  adsorbed  water),  temperature  gradients 
and/or  to  a  microporous  substrate.  The  final  thickness  obtained  in 
the  study  of  a  draining  film  of  octane  on  silicon  by  Sujanani  and 
Wayner  [40]  agreed  with  the  results  of  Ingram.  They  concluded  that 
their  relatively  large  cell  became  contaminated  over  the  long  time 
period  required  for  the  study  of  drainage  and  that  nonisothermal 
zones  were  present  in  the  cell.  Using  SF^,  Kayser  et  al.  [35] 

experimentally  confirmed  that  small  temperature  gradients  have  a 
large  effect  on  the  adsorbed  film  thickness.  Truong  and  Wayner 
found  the  same  effect  in  heat  transfer  studies  [41].  In  a  recent 
study  by  Sujanani  and  Wayner  [42]  on  the  near  equilibrium 
thickness  profile  of  an  extended  meniscus,  the  DLP  theory  and  the 
extreme  sensitivity  of  the  system  to  small  temperature  changes 
were  experimentally  confirmed  and  analyzed.  The  DLP  theory  has 
also  been  confirmed  by  Sabisky  and  Anderson  for  Helium  [31].  The 
DLP  theory  for  alcohols  and  water  has  also  been  confirmed. 
Therefore,  we  find  ample  experimental  evidence  for  "thick" 
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adsorbed  wetting  films  which  are  very  sensitive  to  small 
temperature  differences.  Herein,  we  confirm  this  observation  using 
a  unique  nonequilibrium  system. 

2.1.5  Nonequilibrium  Systems  (Experimental) 

The  isothermal  and  the  evaporating  meniscus  thickness  profiles 
in  the  region  controlled  by  capillarity  have  been  measured  as  a 
function  of  heat  flux  using  interferometry  [e.g.,  15,32,33,39,41].  The 
profile  of  a  wetting  liquid  near  the  transition  region,  where  both 
capillary  and  van  der  Waals  dispersion  forces  are  equally  important, 
has  also  been  measured  using  interferometry  and  ellipsometry 
[28,36,40,42].  In  these  cases,  the  thickness  in  the  capillary  region  was 
determined  as  a  function  of  position  by  measuring  the  locations  of  the 
interference  fringes.  The  thickness  in  the  flat  portion  of  the  adsorbed 
film  was  measured  using  ellipsometry.  In  [42]  additional  physical 
insights  concerning  the  stability  of  the  profile  were  obtained  using 
microcomputer  enhanced  video  microscopy,  wherein  each  pixel  was  a 
photometer,  to  obtain  a  continuous  profile  down  to  the  thickness  of  the 
adsorbed  film  and  then  by  using  ellipsometry  to  confirm  the  thickness 
in  the  flat  region.  The  resulting  equilibrium  and  near  equilibrium 
characteristics  confirmed  past  theoretical  predictions  concerning  the 
thickness  of  the  adsorbed  film  [e.g.,  2,36],  the  width  of  the  transition 
region  [e.g.,  10,11],  and  the  effect  of  heat  transfer  on  the  film  thickness 
profile  [e.g.,  8,9,11].  As  demonstrated  in  this  research,  the  advantages 
of  previously  developed  optical  techniques  to  study  the  details  of 
interfacial  processes  are  significantly  enhanced  when  coupled  with  the 
convenience  of  video  image  processing.  Beaglehole  has  used 
microscopic  imaging  ellipsometry  to  measure  the  profiles  of  the 
precursor  of  spreading  drops  of  a  low  vapor  pressure  fluid  (siloxane) 
on  various  substrates  [38].  Previously,  Bascom  et  al.  measured  the 


profile  of  the  leading  edge  of  a  spreading  drop  using  ellipsometry  [28]. 
Herein,  using  a  new  heat  transfer  cell  we  present  experimental  and 
theoretical  results  that  confirm  the  modeling  represented  by  Eqs.(2.5  - 
2.11)  and  presented  in  Section  4. 
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SECTION  3 


EXPERIMENTAL 

A  cross-sectional  diagram  of  the  circular  capillary  feeder 
system  is  shown  in  Fig.  3.1.  Liquid  flows  from  a  reservoir 
(maintained  at  a  constant  level)  as  a  result  of  a  difference  in 
capillary  pressure,  hydrostatic  head  and  disjoining  pressure  (caused 
by  intermolecular  forces)  and  forms  a  meniscus  at  the  edge  of  the 
thinner  gap  at  position  2.  The  contact  line  region,  area  3,  governs 
the  critical  initial  condition  for  the  stability  and  behavior  of  the 
evaporating  meniscus  and  hence  the  performance  of  the  capillary 
feeder  system.  Liquid  evaporated  from  this  region  is  replenished 
by  fluid  flow  from  the  reservoir.  In  the  thinnest  portion  of  the 
meniscus  represented  by  6o  evaporation  does  not  occur. 

The  first  objective  of  this  study  was  to  determine  the  film 
thickness  profile  in  the  absence  of  external  perturbations.  This 
"isothermal"  profile  gives  the  pressure  field  and  the  interfacial 
properties  (i.e.,  A  by  measuring  So)  and  yields  information  on 
sensitivity  to  evaporation/condensation.  The  capillary  feeder 
system  was  designed  to  provide  a  stable  meniscus  with  the  center 
of  the  silicon  wafer  heated  to  induce  evaporation.  The  change  in 
the  pressure  field  was  manifested  as  a  change  in  shape  of  the  film 
thickness  profile  which  was  measured.  This  data  coupled  with  the 
interfacial  characteristics  and  the  macroscopic  temperature  profile 
allowed  the  evaluation  of  transport  processes.  Using  the  capillary 
feeder  system,  the  amount  of  fluid  evaporated  from  the  entire 
meniscus  was  also  obtained.  Although  the  main  objective  of  this 
study  was  the  evaluation  of  the  theoretical  model  at  low 
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Figure  3.1  Cross  Sectional  View  of  Capillary  Feeder. 
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evaporation  rates,  oscillating  meniscus  profiles  at  higher  power 
inputs  were  also  observed. 

3.1  Heat  Transfer  Cell 

An  aluminum  trapezoidal  cell  chamber  was  fabricated  to 
enclose  the  circular  capillary  feeder  system,  which  allowed  for 
measurable  continuous  liquid  feed  with  simultaneous 
interferometric  and  ellipsometric  measurements  of  film  thicknesses. 
A  photograph  of  the  system  is  presented  in  Figure  3.2.  Two  wide 
circular  grooves  with  different  channel  depths  as  shown  in  Figure 
3.1  were  made  by  a  milled  aluminum  washer  pressing  against  a 
smooth  silicon  wafer  with  a  Teflon  encapsulated  0-ring  in  between. 
The  difference  in  groove  depths  maintained  a  steady  flow  of  liquid 
towards  the  center  of  the  silicon  plate  by  capillary  pumping  to 
replenish  the  amount  evaporated.  The  circular  test  area  and  a 
circular  heat  source,  made  by  baking  a  platinum  compound, 
naturally  eliminated  edge  effects.  Small  (0.3  mm)  chromel-alumel 
thermocouples  (not  shown)  were  attached  1.5  mm  apart  on  the 
bottom  of  the  silicon  substrate  to  measure  the  temperature.  Special 
optical  windows  were  made  at  70°  angles  with  the  base  for 
maximum  ellipsometric  sensitivity.  The  capillary  feeder  was 
connected  to  an  outside  reservoir  (a  graduated  cylinder  which 
contained  pure  test  liquid)  by  means  of  a  very  thin  capillary  tube 
which  supplied  fresh  liquid  by  capillary  suction  to  replenish  the 
liquid  evaporated  in  the  test  area.  It  was  very  important  to 
manually  keep  the  liquid  levels  nearly  the  same  in  the  reservoir 
and  the  feeder.  A  small  change  in  the  hydrostatic  head  could 
overflow  or  empty  the  feeder.  For  the  large  dimensions  used  in  the 
experimental  cell,  we  found  the  system  very  sensitive  to  small 


Figure  3.2  Photograph  of  the  Heat  Transfer  Cell. 
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changes  in  pressure.  The  total  rate  of  evaporation  was  calculated  by 
measuring  the  amount  of  liquid  added  to  the  outside  reservoir. 

3.2  Macroscopic  Heat  Transfer 

Representative  temperature  profiles  for  two  input  power 
settings  with  and  without  fluid  in  the  cell  are  presented  in  Figure 
3.3.  From  the  macroscopic  temperature  profiles,  the  total  rate  of 
heat  loss  (conductive,  convective  and  evaporative),  Ql,  front  the 
whole  extended  meniscus  between  Ri  and  R2  was  calculated  for  the 
dry  and  wet  cases  using; 

QL  =  2iiks^Rl  -  R2 ‘^lR2j  (3.1) 

Here  Ri  is  chosen  as  R  =  10  mm,  beyond  which  the 
evaporation  rate  is  small.  R2  is  chosen  as  R  =  7  mm,  because  for  R  < 

7  mm  only  a  thin  adsorbed  film  exists  where  evaporation  is 
presumed  to  be  negligible.  For  the  dry  case  QL.dry  is  equal  to  Qlc. 
where  QlC  is  only  the  conductive  and  convective  losses.  But  for  the 
wet  (with  liquid)  case  QL,wet  is  related  to  QlC  by  the  following 
relation; 

QL,wet=  Qlc  +  Qev  (3.2) 

Where  Qev  is  the  evaporative  heat  loss.  If  we  assume  that  the 
heat  loss  due  to  conduction  and  convection  mechanisms  does  not 
vary  significantly  between  the  cases  when  the  plate  is  dry  and 
when  the  liquid  is  introduced  to  form  a  meniscus,  the  evaporative 
heat  losses  may  be  calculated  as; 

Qev  =  2*ks(Rl(^lR,  „,„  -^lR,,d,y) 

■^2  idR'*^2,wet'^^R2,dry^^ 
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TEMPERATURE, 


This  evaporative  heat  loss  was  successfully  compared  to  the 
one  calculated  using  the  measured  macroscopic  liquid  input  rate  at 
the  reservoir,  Fr: 

Qt  =  Trhfg  (3.4) 

The  temperatures  obtained  from  the  thermocouples  were 
plotted  against  the  radial  distance,  R,  to  construct  the  temperature 
profiles  for  different  power  inputs  for  both  dry  and  wet  plates.  The 
finite  difference  Levenberg-Marquardt  algorithm  was  used  to  fit  the 
temperature  data  with  radial  distance.  The  functional  relationship 
was  then  differentiated  to  find  dT/dR  at  different  R  values.  Finally, 
using  Eq.  (3.3),  Qev,  the  integral  value  of  the  evaporation  rate  for 
different  power  inputs  was  obtained.  The  evaporative  heat 
required  was  also  determined  from  the  amount  of  liquid  supplied  to 
the  reservoir  using  Eq.  (3.4).  In  Table  1  the  values  of  Qev  obtained 
using  these  different  methods  are  presented.  The  slightly  higher 
value  of  Qev.  predicted  at  the  highest  power  inputs  by  the  model, 
may  be  a  result  of  our  assumption  that  QlC  is  the  same  for  both  dry 
and  wet  cases.  Conceivably,  QLC,wet  is  higher  than  QLC,dry  because 
of  the  presence  of  vapor.  Hence  the  left  hand  side  of  Eq.  (3.3) 
should  be  Qev  +(QLC,wet  -  QLC.dry)- 

Table3.1.  Integral  Evaporation  Rate  as  a  Function  of  the  Heat  Input  Rate 


Power  (watts) 

(watts)  i 

Qev 

Liquid:  Heptane 

1.0 

0.041 

0.042 

2.0 

0.075 

0.060 

3.5 

0.126 

0.120 

4.0 

0.170 

0.190 

Liquid:  R113 

0.25 

0.033 

0.030 

0.35 

0.038 

0.041 

0.50 

0.045 

0.049 

iEq.  (3.4); 

2Eq.  (3.3) 
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We  conclude  that  a  heat  balance  was  obtained  since  Q-p  agrees 
reasonably  well  with  This  confirms  that  substantial 

evaporation  was  present.  Only  a  small  percentage  of  the  input 
power  is  used  to  evaporate  the  liquid  because  of  the  substantial 
values  of  the  convective  and  conductive  terms. 

3.3  Optical  Measurements 

3.3.1  Ellipsometry 

Ellipsometry  uses  the  interference  of  elliptically  polarized 
monochromatic  light  to  analyze  and  measure  the  film  thickness  of 
thin  films.  It  is  a  very  powerful  yet  sensitive  method,  and  is 
applicable  to  a  wide  variety  of  situations.  In  ellipsometry,  the 
change  in  state  of  polarization  due  to  reflection  is  measured  in  terms 
of  properties  of  the  reflecting  surface.  From  a  clean,  film-free 
surface,  the  analysis  of  elliptically  polarized  reflection  can  be  used 
to  determine  the  optical  constants  of  the  surface,  namely  the 
refractive  index  and  the  extinction  coefficient.  For  a  surface  having 
a  thin  film,  given  the  properties  of  the  substrate,  one  can  determine 
the  film  thickness  and  if  the  film  is  thick  enough  (5  >  30  nm),  the 
refractive  index  of  the  thin  film. 

A  schematic  diagram  of  the  experimental  set-up  for  the 
ellipsometer  is  shown  in  Figure  3.4.  A  conventional  null 
ellipsometer  with  the  configuration  of  light  source,  polarizer,  fixed 
compensator,  sample,  analyzer,  microscope  objective  lenses  and  a 
photodetector  was  used  to  determine  the  film  thickness  of  the  flat 
adsorbed  film  which  forms  ahead  of  the  meniscus.  The  light  source 
was  a  Hughes  Model  3022H  He-Ne  laser  of  5  mW  power  which 
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SYMBOLS 


A:  Analyzer 

C  Compensator 
CPL.  Circularly  Polarized  Light 
D  Detector 
F:  Filter 

FCG;  Frame  Code  Generator 
hRV  Image  Resolution  VCR 
I  PM  Image  Processing 
Microcomputer 
L;  Helium-Neon  Laser 


M;  Microscope 
MUL.  Monochromatic 
Unpolanzed  Light 
P:  Polarizer 

PM:  Photomultiplier 
PPL:  Plane  Polarized  Light 
VA:  Variable  Aperture 
SS:  Scanning  Stage 
VC;  CCD  Video  Camera 


Figure  3.4  Schematic  Diagram  of  Experimental  Set-up. 
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emitted  unpolarized,  monochromatic  light  at  a  fixed  wavelength  of 
6328A  When  plane  polarized  light  is  incident  on  a  solid  surface,  the 
reflects  i  light  becomes  elliptically  polarized.  The  changes  in 
polarization,  i.e.,  the  relative  amplitude  (y)  and  phase  angle  (6), 
depend  on  the  optical  constants  of  the  solid,  the  angle  of  incidence  and 
the  wavelength  of  light.  If  there  is  a  thin  film  adsorbed  on  the 
surface,  H'  and  6  will  also  depend  on  the  optical  constants  and 
thickness  of  the  film.  A  computer  program  based  on  the  algorithm  of 
McCrackin  was  used  to  calculate  the  film  thickness  from  '¥  and  6  which 
were  measured  in  two  optical  zones  and  zone  averaged  [36,51].  The 
measuring  precision  of  our  ellipsometry  system  was  checked  prior  to 
the  experiment  by  measuring  the  optical  constants  (refractive  index) 
of  a  bare  silicon  wafer.  (Measured  optical  constants  for  silicon  ns=3.85 
±  0.01,  ks=0.05  ±  0.02;  literature  values  ns=3.882,  ks=0.019  [36].) 


3.3.2  Interferometry 


The  interference  phenomena  associated  with  monochromatic 
light  reflected  from  a  tapered  film  were  used  to  determine  the 
profile  of  the  capillary  meniscus  in  the  thickness  range  5>  0.1  pm. 
Alternating  constructive  and  destructive  interference  occur  when 
coherent  light  is  reflected  from  the  two  interfaces  of  the  thin  liquid 
film  and  recombine.  Destructive  interference  occurs  when  the  optical 
paths  of  the  two  reflected  beams  are  such  that  they  are  out  of  phase 
by  7t.  In  our  case  ny  <  ni  <  ns;  therefore,  there  is  a  phase  shift  equal 
to  It  at  each  interface.  The  relevant  equation  for  the  film  thickness 
associated  with  the  minima  is: 


5 


(2L  +  1) 
4n/ 


L  =  0,1,2,... 


(3.5) 


where  n  is  the  refractive  index,  subscripts  1,  v  and  s  denote  liquid, 
vapor  and  solid  respectively,  \  is  the  wavelength  of  the 
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monochromatic  light  and  L  is  the  order  of  the  destructive  (dark) 
fringes.  We  used  a  high  power,  long  working  distance  microscope 
(Wild  Model  M420)  to  view  the  naturally  occurring  interference 
fringes  in  the  thicker  portion  of  the  thin  film.  A  picture  showing 
representative  sample  is  presented  in  Figure  3.S.  Coherent  light  (X  = 
S2S  nm)  was  used  at  normal  incidence  and  the  reflected 
interference  patterns  were  captured  using  a  CCD  camera  (Sony 
HVM200). 

3.4  Image  Processing 

To  automate  data  acquisition,  improve  data  resolution  and 
enhance  data  analysis,  we  analyzed  the  interferometric  images 
using  an  image  processing  system.  The  first  step  in  digital  image 
processing  is  the  conversion  of  the  optical  image  into  a  form  that 
can  be  stored  in  computer  memory.  This  conversion  is  performed 
by  a  light  sensitive  system  known  as  an  optical  digitizer,  which 
produces  coded  numbers  that  are  a  measure  of  light  intensity  in  the 
optical  image.  This  process  is  called  digitizing,  and  the  stored 
numerical  representation  of  the  original  optical  image  is  called  a 
digital  image.  The  digital  image  can  be  manipulated  like  any  other 
data  stored  in  computer  memory.  All  such  manipulations  can  be 
placed  in  one  of  two  general  categories;  image  analysis  or  image 
enhancement.  The  term  image  analysis  refers  to  the  use  of  digital 
computers  to  derive  numerical  information  regarding  selected 
image  features,  such  as  lengths,  areas,  intensity  distribution,  etc. 
Image  enhancement,  on  the  other  hand,  is  a  procedure  for 
manipulating  the  video  signal  to  sharpen  or  otherwise  improve  the 
video  picture.  Routines  for  image  enhancement  differ  from  those  of 
image  analysis  in  that  the  object  of  the  processing  is  not  to  produce 
descriptive  information  about  the  image,  but  rather  to  manipulate 
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the  image  in  such  a  way  that  it  is  more  useful  to  a  human  observer. 
We  used  both  the  techniques  in  our  work. 

Briefly,  the  system  consisted  of  a  16  MHz  IBM'PC286  clone 
computer  equipped  with  real  time  frame  grabber  and  frame 
processor  boards  (DT-2851  and  DT-2858  from  Data  Translation 
Inc.).  The  image  captured  from  the  microscope  through  a  CCD 
camera  was  digitized  into  490  (vertical)  x  S12  (horizontal)  space 
pixels  and  assigned  one  of  2S6  possible  grey  values  representing  an 
intensity  from  0  (black)  to  25S  (white).  Thus  the  image  area, 
already  magnified  considerably  by  the  microscope,  was  divided  into 
microscopic  pixels  which  acted  as  individual,  simultaneous  light 
sensors.  Each  picture  provided  data  on  one  frame;  i.e.,  l/30th  of  a 
second  in  real  time.  The  interferometric  fringes  formed  were 
parallel  because  the  radial  distance  from  the  center  [~8.8  mm]  was 
large  compared  to  the  meniscus  width  [~0.2  mm].  When  viewed  at 
high  magnification,  (area  viewed  -0.3  mm  x  0.3  mm]  the  fringes 
were  almost  invariant  in  a  direction  perpendicular  to  the  radius  and 
could  therefore  be  averaged  spatially  to  improve  the  signal/noise 
ratio.  The  light  level  was  optimized  to  prevent  pixel  saturation 
(excessive  light)  and  to  ensure  linearity  of  the  analog  to  digital 
conversion.  Light  below  a  certain  level  drastically  increased  the 
signal/noise  ratio.  For  our  system  the  best  combination  of 
resolution,  light  level,  signal/noise  ratio  and  magnification  was 
found  by  adding  a  4.0  X  effective  magnification  in  the  tube.  Each 
pixel  represented  a  diameter  of  0.62S  pm.  Transient  data  were 
captured  by  continuously  recording  on  a  high  quality  U-Matic  Sony 
V0'9600  VCR.  The  VCR  had  a  B/W  resolution  of  more  than  300 
lines  and  was  connected  through  a  Sony  FCG-700  Frame  Counter, 


which  superimposed  each  frame  with  a  permanent  identification 
number. 


3.5  System  Operation 

Silicon  was  used  as  the  substrate  because  of  its  reflective, 
polished,  smooth,  easily  available  and  reproducible  surface.  We 
note  that  a  native  oxide  about  3~S  nm  thick  is  always  present  on 
the  surface.  The  exact  composition  of  this  layer  depends  on  the 
prior  history  of  the  wafer.  All  cleaning  procedures  were  done  inside 
a  class- 100  clean  hood.  The  cell  was  first  immersed  in  pure  ethanol 
(absolute-200  proof)  for  6  to  7  hours.  The  cell  pans  were  again 
rinsed  with  fresh  alcohol  and  then  dried  in  an  oven  at  180°C. 

Extreme  care  was  taken  to  ensure  that  no  dust  particles  or 
impurities  could  come  in  contact  with  the  silicon  wafer.  After 
drying,  the  cell  parts  were  rinsed  with  pure  test  fluid  and  blow 
dried  with  extra-dry  grade  nitrogen.  The  assembled  cell  was 
mounted  on  the  stage  of  a  high  power  long  working  distance  'WILD 
microscope.  The  system  was  left  to  equilibrate  with  the 
surroundings  for  about  3  hours  before  taking  data.  An  optical 
interference  pattern  representing  the  thickness  profile  was  readily 
observed:  an  example  is  given  in  Figure  3.S.  The  pattern  was  then 
analyzed  to  determine  the  film  thickness  profile  as  detailed  in  the 
next  section.  The  temperature  distribution  along  the  wafer  and  the 
macroscopic  liquid  flow  rate  were  also  noted. 

3.6  Image  Analysis 

From  a  picture  a  plot  of  the  pixel  grey  value  versus  distance 
was  extracted.  The  grey  value  at  each  pixel  is  a  measure  of  the 
reflectivity  which  changes  with  film  thickness.  In  particular  for  a 
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thin  liquid  film  of  refractive  index,  ni,  on  a  solid  surface  of 
refractive  index,  ns. 


RL  (Reflecrivity)  = 


g  +  pcos26i 
K  +  pcosiei  ’ 


01  =  2jeni8A;  o  =  ;  3  =  2rir2  ; 


(3.6) 


K 


ni  -  ny  ns  -  ni 

ni  +  nv  '  ns  +  oi 


where  nv  is  the  refractive  index  of  vapor  and  X  is  the  wavelength  of 
light  in  vacuum.  As  is  evident  from  Eq.  (3.6),  the  reflectivity 
undergoes  a  cyclic  change  with  increase  in  film  thickness. 


RLmax  = 

+  T2  Y 

,1  +  rir2j 

a  +  3 
“  K  +  3  ’ 

01  =  0,Lii 

(3.7) 

RLmin  = 

ft\  -  r2  ^ 
VI  -  rir2j 

0-3 
"k  -  3  ’ 

n  1 

01  =  (L  +  2^^* 

(3.?) 

The  computer  program  scanned  the  peaks  and  valleys  and 
filtered  the  noise  from  the  real  peaks/valleys.  It  then  interpolated 
peak/valley  envelopes  and  by  analyzing  the  relative  reflectivity  of 
any  given  pixel  with  respect  to  these  (dark  and  light  pixel 
envelopes),  determined  a  film  thickness  at  every  pixel.  Thus, 


Gjx)  -  Gmin(x) 

*  Gmax(x)  -  Gmin(x) 


(3.9) 


201  =  Cos-1 


(&  +  K 


1  •  2  G(x))l 


{2G(x)  -  1)  -  K  j 


(3.10) 


where  Gmin(x)  and  Gmax(x)  are  the  interpolatory  envelopes  to  the 
various  order  minima  and  maxima  (except  for  the  first  dark  fringe. 
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where  linear  extrapolation  and  matching  with  ellipsometric  data 
were  done).  By  knowing  0i  and  the  order  'L'  of  the  fringe,  6  was 
found. 

3.7  Experimental  Results 

Two  examples  of  the  measured  thickness  profiles  are  shown  in 
Fig.  3.6.  The  profiles  were  shifted  on  the  x  axis  to  the  same  relative 
position  at  a  fixed  thickness.  The  adsorbed  film  thicknesses  in  the 
flat  region  were  So  =  7.5  nm  and  So  =  15  nm.  The  profile  for  So  = 

15nm  can  be  termed  as  a  "very  near  equilibrium  case  with  very 
small  evaporation, "  whereas  So  =  7.5  nm  is  a  "near  equilibrium  case 

with  small  evaporation."  Subsequent  analysis  and  discussion  in  this 
report  will  clarify  these  terminologies.  Figure  3.7  shows  thickness 
profiles  obtained  during  evaporation  with  controlled  rates  of  energy 
input,  Q.  These  two  figures  clearly  show  that  the  isothermal  profiles 
are  more  spread  out  and  that  the  thickness  of  the  adsorbed  film,  6o, 

decreases  with  an  increase  in  power  input  (surface  temperature). 

In  an  isothermal  horizontal  system  of  spreading  liquid  on  a  solid 
substrate  at  the  exit  of  a  capillary  the  curvature  should  remain 
constant  in  the  region  where  dispersion  forces  can  be  neglected.  So 
the  film  profile  in  this  range  approximates  a  parabola  and  a  plot  of 
51/2  vs  X  is  a  straight  line.  If  the  curvature  is  constant,  it  is  related 
to  the  slope  of  this  line  as  Kc  =  2(d6l/2/dx)2.  The  experimental 
values  of  51/2  ^re  also  given  in  Figs.  3.6  and  3.8.  For  the  case 
where  Sq  -  15  nm,  51/2  vs  x  is  nearly  a  straight  line  which  means 
the  system  was  very  close  to  isothermality,  but  for  5©  =  7.5  nm  the 
line  bends  downward,  showing  the  presence  of  a  curvature  gradient 
and  hence  some  evaporation.  Although  the  heater  was  off,  Q  =  0,  a 
small  amount  of  heat  conduction  from  the  surroundings  caused  the 
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Figure  3.6  Film  Thickness  Profile  for  Heptane  on  Silicon 
with  a  Native  Surface  Oxide.  (Very  near 
equilibrium  with  very  small  evaporation. 

Q  »  O.OW,  5o  »  15  nm  and  near  equilibrium 
with  small  evaporation,  Q  »  O.OW.  8o  »  7.5  nm). 
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Figure  3.7  Film  Thickness  Profiles  for  Heptane  on  Silicon 
with  a  Native  Surface  Oxide.  (Very  near 
equilibrium  with  very  small  evaporation. 

5<j  =  15  nm  and  evaporation  with  controlled 
rate  of  heat  inputs.  5o  *  6.2  nm  and  5.6  nm). 
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Figure  3.8  vs.  X. 
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phase  change.  Other  film  thickness  profiles  are  presented  in  Figs 
3.9-3.22. 

To  compute  the  meniscus  curvature,  the  profile  data  need  to 
be  differentiated  twice,  but  magnifications  of  experimental  errors 
need  to  be  controlled  at  every  differentiation  step.  In  References 
[42,32],  the  process  of  cubic  spline  fitting  and  numerical 
differentiation  was  used.  Herein,  we  fit  the  data  with  a  theoretical 
model  as  described  in  Section  4.  The  current  procedure  seems  to  be 
more  accurate  because  of  the  location  of  the  maximum  curvature 
and  the  effect  of  surface  tension.  Further,  the  comparison  between 
the  theoretical  and  experimental  thickness  profiles  is  easy  io 
evaluate. 


37 


^0  so  100  ISO  200  2S0 


RELATIVE  DISTANCE,  X. 

Figure  3.9  Film  Thickness  Profiles  of  Heptane  on  Silicon 
with  a  Native  Surface  Oxide.  (Very  near 
equilibrium  with  very  small  evaporation, 

5o  =  15  nm  and  evaporation  with  controlled 
rate  of  heat  inputs,  5o  =  4.8  nm  and  4.4  nm). 


Figure  3.10  vs.  X. 
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Figure  3.11  Film  Thickness  Profiles  of  R113  on  Silicon 
with  a  Native  Surface  Oxide.  (With  no  heat 
input.  So  -  11.5  nm  and  evaporation  with 
controlled  rate  of  heat  inputs.  So  =  4.8  nm 
and  4.0  nm). 
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Figure  3.12  5>/2  vs.  X. 
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Figure 
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.13  Film  Thickness  Profiles  of  R113  on  Silicon 
with  a  Native  Surface  Oxide.  (Evaporation 
with  controlled  rate  of  heat  inputs.  So  = 

4.8  nm,  4.S  nm  and  4.0  nm). 
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Figure  3.14  vs.  X. 
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Figure  3. IS  Film  Thickness  Profiles  of  Propanol  on 

Silicon  with  a  Native  Surface  Oxide.  (With  no 
heat  input.  So  ^  13.6  nm  and  evaporation 
with  controlled  rate  of  heat  inputs, 

So  -  3.6  nm  and  2.9  nm). 
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Figure  3.16  Si/2  vs.  X. 
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Figure  3.17  Film  Thickness  Profiles  of  Propanol  on 
Silicon  with  a  Native  Surface  Oxide. 
(Evaporation  with  controlled  rate  of  heat 
inputs,  6o  =  6.5  nm,  5.2  nm  and  4.7  nm). 
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Figure  3.18  51/2  vs.  X 
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Figure  3.19  Film  Thickness  Profiles  of  Butanol  on 

Silicon  with  a  Native  Surface  Oxide.  (With 
no  heat  input.  So  =  17.2  nm  and  evaporation 
with  controlled  rate  of  heat  inputs. 

So  =  5.9  nm  and  4.9  nm). 


^0  50  too  150  200  250 

RELATIVE  DISTANCE,  X,  fim 


Figure  3.20  5i/2  vs.  X. 
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®0  50  100  150  200  250 

ftCUTIVE  distance,  X, 

Figure  3.21  Film  Thickness  Profiles  of  Butanol  on 

Silicon  with  a  Native  Surface  Oxide.  (With 
no  heat  input,  5  o  ^  17.2  nm  and  evaporation 
with  controlled  rate  of  heat  inputs, 

5o  =  5.4  nm  and  4.5  nm). 


REUTIVE  distance,  X,  pim 
Figure  3.22  vs.  X. 
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SECTION  4 


THEORETICAL 

The  immediate  theoretical  precursors  to  this  work  are  the 
studies  of  Wayner,  et  al.  [9],  Moosman  and  Homsy  [11],  and  Wayner 
[19].  A  substantial  portion  of  the  material  in  the  initial  theoretical 
paper  associated  with  this  project  [19],  which  developed  the 
connection  between  the  film  tension  and  the  following  equations, 
will  not  be  included  herein.  It  was  included  in  the  interim  repon 
(WRDC-TR-90-2008).  The  current  work  includes  all  of  the  effects 
considered  by  Wayner,  et  al.  [9,19]  and  Moosman  and  Homsy  [11] 
but  the  nonlinearity  of  the  combined  effects  is  dealt  with 
numerically.  Consider  a  stationary  thin  film  at  the  leading  edge  of  a 
thicker  pool  of  a  pure  liquid  (also  stationary)  on  a  horizontal 
substrate.  The  liquid  wets  the  solid  substrate  and  therefore  at 
steady  state  the  thin  Him  extends  a  large  distance  from  the  bulk 
pool.  This  configuration  is  shown  in  Figure  4.1.  Because  the  film  is 
thin  it  is  strongly  influenced  by  capillary  effects  and  by  adsorption 
forces.  In  this  study  the  adsorption  forces  are  restricted  to  the 
London-van  der  Waals  dispersion  force  and  the  film  is  assumed  to 
be  flat  enough  that  the  hydrodynamics  are  governed  by  lubrication 
phenomena.  Thus  the  pressure  in  the  liquid  does  not  vary  in  the 
direction  perpendicular  to  the  substrate. 

Liquid  flows  in  the  steady  film  because  of  evaporation 
resulting  from  superheating.  The  local  evaporation  rate  is  affected 
by  the  local  superheat  and  the  local  film  pressure.  The  film 
pressure  reflects  the  adhesion  of  the  liquid  to  the  solid  substrate 
and  the  capillary  pressure.  The  tendency  of  the  liquid  to  evaporate 
is  reduced  by  adhesion  and  a  negative  capillary  pressure.  (A 
negative  capillary  pressure  is  caused  by  a  thin  film  which  is  concave 


Figure  4.1 


Theoreiical  Configuration. 
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towards  the  vapor  phase.)  The  local  superheat  is  the  temperature 
difference  between  the  liquid-vapor  interface  and  the  vapor  phase. 
The  vapor  phase  is  assumed  to  be  isothermal  and  pure. 

Heat  flows  from  the  substrate  through  the  film  to  the  liquid- 
vapor  interface  where  evaporation  occurs.  Since  the  substrate  is 
assumed  to  be  isothermal,  and  the  film  has  a  resistance  to  heat 
conduction,  we  expect  that  the  liquid-vapor  interface  of  the  thicker 
portions  of  the  film  will  be  less  superheated  than  the  thinner 
portions.We  assume  that  the  film  extends  an  infinite  distance  from 
the  pool.  If  the  film  is  thin  enough,  the  liquid-solid  intermolecular 
force  field  is  sufficiently  strong  to  keep  the  film  from  evaporating 
even  though  it  is  superheated.  Therefore,  we  expect  the  evaporation 
rate  to  be  highest  where  the  film  thickness  is  moderate.  Far  from 
the  liquid  pool  the  film  does  not  evaporate  at  all.  The  variation  of 
surface  tension  with  temperature  is  neglected. 


The  literature  provides  relationships  which  are  the  building 
blocks  of  the  model  of  the  horizontal  evaporating  thin  film  [e.g., 
9,11,19].  The  film  thickness  is  deHned  to  be  5,  and  the  coordinate 
along  the  solid  substrate  is  x  as  shown  in  Figure  (4.1).  The  liquid 
pressure,  P/,  is  related  to  the  vapor  phase  pressure,  Py.  by; 


Pi-Pv 


A  d^5 
“  53  ■  ®  dx2 


(4.1) 


where  6xA  is  the  Hamaker  constant  which  is  negative  for  a 
spreading  liquid,  and  o  is  the  surface  tension.  The  first  term  in  Eq. 
(4.1)  is  called  the  disjoining  pressure  and  represents  the  body  force 
on  the  liquid  due  to  van  der  Waals  dispersion  phenomena. 

According  to  Eq.  (4.1),  the  adhesion  effect  is  larger  for  a  thinner  film. 
The  expression  used  for  disjoining  pressure  is  for  a  flat  film,  so  we 
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require  the  film  to  be  fairly  flat.  According  to  a  two-dimensional- 
force  balance  for  a  nonevaporating  static  film  by  Kralchevsky  and 
Ivanov  [53],  the  surface  tension  in  a  thin  film  is  not  constant. 
Applying  their  analysis  to  a  fairly  flat  film  we  find  the  variation  of 
the  surface  tension  to  be  inconsequential.  The  simplified  form  of 
curvature  is  used  because  the  film  is  assumed  to  be  fairly  flat.  The 
precise  curvature  is: 


K  = 


d28 
_ dx^ 

[1 


(4.2) 


therefore,  the  approximate  form  used  is  only  valid  if: 


«  1 


(4.3) 


Lubrication  theory,  which  is  valid  only  if  Equation  (4.3)  is  satisfied, 
relates  the  mass  flow  rate  in  the  thin  film  to  the  pressure  gradient 
in  the  direction  of  flow: 

where  r(x)  is  the  mass  flow  rate  and  v  is  the  kinematic  viscosity. 
Equations  (4.1)  and  (4.4)  show  that  a  given  thin  film  shape,  6(x), 
gives  rise  to  a  given  mass  flow  rate,  r(x). 

Following  Wayner  [19],  the  evaporative  flux  is  modelled 
according  to  the  expression: 

m  =  a(T/v  -  Tv)  +  b(P/  -  Pv)  (4.5) 

where 

m  is  the  evaporative  flux,  defined  in  the  direction  away  from 
the  film: 
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T/v  is  the  temperature  of  the  liquid-vapor  interface,  and 


Tv  is  the  temperature  of  the  vapor: 


a  and  b  are  functions  of  the  system  temperature  and  physical 
properties  and  are  defined  as: 


b  =  C 


f  M 

^  vMAhm 

|,2itRT/vJ 

1  RTvT/v 

r  M  y/2 

fV/Pv) 

^2ttRT/v  J 

[rT/v  j 

(4.6) 

(4.7) 


where 

C  is  the  accommodation  coefficient  taken  to  be  2.0; 

M  is  the  molecular  weight  of  the  evaporating  material; 

R  is  the  universal  gas  constant; 

Pv  is  the  bulk  vapor  pressure  of  the  evaporating  material  (e.g., 
uninfluenced  by  disjoining  pressure,  surface  tension,  etc.)  at 
temperature  Tv; 

V/  is  the  molar  volume  of  the  bulk  material  at  temperature 
T/v; 


Ahfn  is  the  enthalpy  of  vaporization  of  the  material  per  unit 
mass,  at  temperature  T/v. 

In  practice  the  temperature  T/v  is  unknown;  however,  if  the 
temperature  differences  between  the  substrate,  vapor,  and  interface 
are  small  compared  to  the  absolute  temperature  of  the  substrate, 
then  the  known  substrate  temperature,  Ts,  may  be  substituted  for 
the  unknown  T/v.  Equation  (4.S)  says  that  evaporation  is  promoted 
by  a  superheat  and  hindered  by  low  liquid  film  pressures.  This 
model  implies  that  the  vapor  phase  does  not  present  a  significant 
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resistance  to  evaporation.  The  temperature  of  the  liquid-vapor 
interface  is  related  to  the  temperature  of  the  substrate  through  the 
one-dimensional-conduction  heat  transfer  solution  for  the  film,  valid 
if  Equation  (4.3)  is  satisfied. 


g  (Ts  -  T/v)  =  m  Ahm 

where  k  is  the  thermal  conductivity  of  the  liquid  in  the  film. 
Following  Moosman  and  Homsy  [11],  Equation  (4.8)  may  be 
combined  with  (4.5)  to  eliminate  T/v  in  favor  of  Tj: 


(4.8) 


nt  -  *  u  [s(Ts  ■  Tv)  +  b(P/  -  Pv)] 

W  -1^6 


(4.9) 


This  equation  clearly  demonstrates  the  direct  effect  of  film  thickness 
on  the  evaporation  rate.  A  recent  numerical  study  by  Stephan  and 
Busse  also  demonstrated  the  importance  of  conductive  resistance  [58]. 

The  local  evaporation  rate  is  linked  to  the  flow  rate  in  the  film 
through  a  material  balance: 


=  -  m 


(4.10) 


The  coupled  differential  equations  are  thus  Eq.  (4.1)  and: 


J_4_r  3  dj^ 

3v  dx  1  dx 


dPi]  _ 


aAhm  . 
1  —5 


[aAT  b(P/-Pv)] 


(4.11) 


where: 


AT  =  Ts  -  Tv 


(4.12) 


The  variables  are  nondimensionalized  in  the  following  way. 
For  reference,  consider  a  strictly  flat  film: 


5(x,t)  =  6o 


(4.13) 
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There  is  no  pressure  gradient  in  such  a  film  according  to  Equation 
(4.1).  Therefore  the  right-hand  side  of  Equation  (4.11)  must  be  zero. 
Hence, 


b(P/  •  Pv)  =  -HAT 


(4.14) 


Combining  this  result  with  Equations  (4.1)  and  (4.13)  we  find: 


A  a  « 

^3  =  -  h  iT 


(4.15) 


If  So  is  so  chosen,  the  flat  film  is  a  solution  to  the  governing 
Equations  (4.1,  4.11).  The  reference  thickness  is  thus  fixed  and  the 
reference  pressure  is  the  magnitude  of  the  pressure  of  such  a  film. 
This  pressure  is  given  by  (4.14)  as: 


°  =  w 


(4.16) 


The  dimensionless  pressure  is  defined  as: 
^  =  (P/-Pv)/no 


(4.17) 


The  scaling  of  x  is  now  stipulated  by  Equation  (4.11)  as: 


-A 

vaAT 


(4.18) 


(Note:  A  is  negative  for  a  spreading  film  and  m'4  =  a  (T^  -  Ty)  = 


aAT,  the  ideal  evaporative  mass  flux.) 


The  dimensionless  position  and  film  thickness  are: 


%  S  X// 


(4.19) 


q  a  S/So 


(4.20) 


The  dimensionless  form  of  Equations  (4.11)  and  (4.1)  is: 

1 


i  i  (,  di', 


3  d? 


1  +  KTl 


(1  +4>) 


(4.21) 


and: 

1  d2q 

'  1)3  ■  ^  d42 

where  the  parameters  k  and  e  are: 
aAhmSo 


and: 

o5obv 

e  =  — : — 

-A 


(4.22) 


(4.23) 


(4.24) 


The  parameter  x  is  a  measure  of  the  importance  of  the  resistance  of 
the  film  to  thermal  conduction.  The  parameter  e  is  a  measure  of  the 
importance  of  capillary  pressure  effects  relative  to  disjoining 
pressure  effects. 

We  impose  two  far  field  conditions  .  The  first  condition  is  for 

Ti  1  (4.25) 


(|)  ->  -1 


(4.26) 


The  film  thickness  is  assumed  to  asymptotically  approach  the  non¬ 
evaporating  thin  film  thickness  mentioned  previously.  The  second 
condition  is  for  ^  -*  -««: 


n  -♦  oo  _  (4.27) 

0  “0m  (4.28) 
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NUMERICAL  SOLUTION 


The  governing  equations  (4.21,  4.22)  were  transformed  with 
the  substitution: 

X  =  (4.29) 

Ve 


to  yield: 

3  dX  dX  J  ^  ®  1  +  K  n 


(4.30) 


and: 


<D  =  . 


n3  dX2 


(4.31) 


The  far-field  behavior  described  by  Equations  (4.2S)  and 
(4.26)  may  be  further  described  by  analysis  of  Equations  (4.30)  and 
(4.31)  for  r|  near  1  and  4  near  -1.  The  result  is: 


Ti  1 


+  po  exg  (■W  7^  ]*  !*>  «p(-^''3) 


(4.32) 


and: 


•1  +3po  1 


exp  -X 


V  K  +  1 


(4.33) 


as  X  — » 


This  result  provides  four  initial  conditions  to  the  Adams  method 
NAG  routine  IX)2CBF.  As  written  these  are  more  like  "final 
conditions"  than  initial  conditions;  however,  the  independent 
variable  in  the  code  is  actually  -X  and  so  Equations  (4.32)  and  (4.33) 
provide  initial  conditions.  A  value  of  -X  was  chosen  (typically  -2.8), 
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dll  d0 

r|,  -  ♦  and  -  ~  were  calculated,  and  the  routine  was  called. 


The 


value  of  -X  must  be  sufficiently  negative  so  that  r|  is  close  enough  to 

1  and  to  -1  that  the  expressions  (4.32)  and  (4.33)  are  a  valid 

dri  d^ 

approximation  of  n  and  <».  The  analytical  expressions  for  ~  and  ~ 


are  simply  the  derivatives  of  ti  and  0.  Therefore,  the  solutions  rj  and 
(|>  are  functions  of  e,  k,  3o.  ^d  |3  i. 


The  numerical  solution  was  compared  to  the  experimental 
results  in  the  following  way.  The  flat  ponion  of  the  film,  Sq,  was 
measured  with  Ellipsometry.  Using  the  measurement  of  So,  k  was 
calculated  with  Equation  (4.23).  A  value  of  e  was  selected  and  then 
the  ratio  of  3o  to  Pi  was  varied  while  one  of  them  (typically  Pi)  was 
held  constant.  This  ratio  Po/Pi  was  varied  until  the  two  end  points 
of  the  data  were  matched  by  the  numerical  result.  If  the  value  of  e 
was  too  large  the  numerical  solution  curved  too  much  relative  to  the 
data  (the  points  between  the  end  points),  lagging  below.  If  the 
value  of  £  was  too  small  the  numerical  solution  was  too  straight  and 
the  data  between  the  end  points  lagged  below  the  solution.  In  this 
way  the  "best"  value  of  e  was  selected.  The  reason  Pi  was  held 
constant  is  that  the  solutions  r\  and  $  are  invariant  with  respect  to 
translation  of  X,  as  shown  by  Equations  (4.30)  and  (4.31).  That  is  if 
^(X)  satisfies  (4.30,  4.31)  so  does  $(X  +  c)  where  c  is  any  constant. 

The  computer  program  is  given  in  Appendix  C.  The  comparison 
between  data  and  model  was  performed  on  a  dimensional  basis  and 
is  pictured  in  Figures  4.2  to  4.24.  In  these  results,  Q  is  the  power 
input  setting,  So  is  the  measured  value,  k  depends  on  5o,  and  e  is  the 
value  of  the  unknown  being  determined. 
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FILM  THICKNESS. 


RELATIVE  DISTANCE,  X,  /zm 


Figure  4.2  Comparison  of  Theory  and  Experiments 

(Q  =  O.OW,  5o  =  15  nm.  K  =  6.48  x  10-2,  e  =  2.5, 
Po  =  -1.800  X  10-4,  p,  =  1..0  X  10-2). 
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FILM 


20  40  60  80 

RELATIVE  DISTANCE,  X,  /um 


Figure  4.3  Comparison  of  Theory  and  Experiments 

(Q  =  O.IW,  6o  =  6.4  nm,  ic  =  2.771  x  10-2,  e  =  1.1, 
Po  =  -6.930  X  10*3,  Pi  =  1.0  X  10*2). 
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FILM  THICKNESS. 


RELATIVE  DISTANCE,  X.  ium 


Figure  4.4  Comparison  of  Theory  and  Experiments 

(Q  =  0.2W.  8o  =  6.2  nm,  k  =  2.678  x  10-2,  e  =  1.5, 
Po  =  -1.712  X  10-3,  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS, 


RELATIVE  DISTANCE.  X.  /urn 


Figure  4.5  Comparison  of  Theory  and  Experiments 

(Q  =  0.3W,  5o  =  5.9  nm,  k  =  2.554  x  10-2,  e  =  1.1, 
Po  =  -6.858  X  10-3,  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS,  6 


RELATIVE  DISTANCE,  X,  yum 


Figure  4.6  Comparison  of  Theory  and  Experiments 

(Q=  0.5W,  5o  =  5.65  nm,  k  =  2.440  x  10-2,  e  =  1.05, 
Po  =  *8.645  X  10-3,  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS 


RELATIVE  DISTANCE.  X,  nrr\ 


Figure  4.8  Comparison  of  Theory  and  Experiments 

(Q  *  O.OW,  6o  =  11.5  nm,  k  =  2.86  x  lO'i,  e  =  6.0, 
Po  =  -6.900  X  10-6,  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS, 


RELATIVE  DISTANCE,  X,  /xm 


Figure  4.9  Comparison  of  Theory  and  Experiments 

(Q  =  0.25W,  8o  =  4.8  nm,  k  =  1.194  x  10*1,  e  =  5.5, 
Po  =  -4.500  X  10-6,  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS. 


RELATIVE  DISTANCE,  X,  /xm 


Figure  4.10  Comparison  of  Theory  and  Experiments 

(Q  =  0.35W,  6o  =  4.5  nm.  k  =  1.193  x  10- 1.  e  =  8.0, 
Po  =  -3.950  X  10-7,  Pi  =  1.0  X  10-2). 
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FILM 


Figure  4.11  Comparison  of  Theory  and  Experiments 

(Q  =  0.5W.  5o  =  4.0  nm,  k  =  1.12  x  lO-l,  e  =  4.0, 
Po  =  -2.630  X  10-5,  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS, 


RELATIVE  DISTANCE,  X, 


Figure  4.12  Comparison  of  Theory  and  Experiments 

(Q  =  O.OW,  8o  =  13.6  nm,  k  =  4.74  x  10*2,  e  =  5.5, 
3o  =  '^-500  X  10-2,  Pi  =  1.0  X  10-2). 
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FILM  T 


RELATIVE  DISTANCE,.  X,  /xm 


Figure  4.13  Comparison  of  Theory  and  Experiments 

(Q  =  0.05W,  6o  =  6.5  nm,  k  =  2.26  x  10-2,  e  =  5.5 
3o  =  -2.100  X  10-6,  p,  =  1.0  X  10-2). 
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FILM  THICKNESS 


RELATIVE  DISTANCE,.  X,  /xm 


Figure  4.14  Comparison  of  Theory  and  Experiments 

(Q  =  0.075W,  6o  =  5.2  nm,  k  =  1.81  x  10*2,  t  =  5.5 
(Jo  =  -2.100  X  10-6,  Pi  =  1.0  X  10-2). 
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FILM  T 


RELATIVE  DISTANCE,  X,  /um 


Figure  4.15  Comparison  of  Theory  and  Experiments 

(Q  =  O.IW,  8o  =  4.7  nm,  k  =  1.64  x  10-2,  £  =  5.5) 
(Jo  =  -2.100  X  10-6,  p,  =  1.0  X  10-2). 
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FILM  THICKNESS, 


RELATIVE  DISTANCE,  X,  txm 


Figure  4.16  Comparison  of  Theory  and  Experiments 

(Q  =  0.1 25W,  5o  =  4.4  nm,  k  =  1.53  x  10*2.  e  =  5.5 
Po  =  -2.100  X  10-6,  pi  =  1.0  X  10-2). 
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FILM  THICKNESS, 


RELATIVE  DISTANCE,  X,  fim 


Figure  4.17  Comparison  of  Theory  and  Experiments 

(Q  =  0.175W,  6o  =  3.6  nm.  k  =  1.44  x  10-2,  e  =  5.5 
Po  =  -2.100  X  lO-fi,  Pi  =  1.0  X  10-2). 
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RELATIVE  DISTANCE,  X, 


Figure  4.18  Comparison  of  Theory  and  Experiments 

(Q  =  0.54W,  5o  =  2.9  nm,  k  =  1.36  x  10-2,  e  =  5.0 
3o  =  -3.825  X  I0-<>.  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS 


RELATIVE  DISTANCE,  X, 


Figure  4.19  Comparison  of  Theory  and  Experiments 

(Q  =  0.77W,  5o  =  2.7  nm,  k  =  1.41  x  10-2,  e  =  i.g 
00  =  '7.285  X  10-4,  0,  =  10  X  10-2). 
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RELATIVE  DISTANCE,  X.  (xm 


Figure  4.20  Comparison  of  Theory  and  Experiments 

(Q  =  O.OW,  5o  =  17.2  nm,  k  =  1.652  x  10*2.  e  =  ll.O 
3o  =  1.000  X  10-7.  3i  »  1.0  X  10-2). 
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FILM  THICKNESS 


RELATIVE  DISTANCE,  X,  /xm 


Figure  4.21  Comparison  of  Theory  and  Experiments 

(Q  =  0.05W,  8o  =  5.9  nm,  k  =  6.063  x  lO'^,  e  =  5.0, 
00  =  -3.300  X  10-6,  p,  =  10  X  10-2). 
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FILM  THICKNESS, 


RELATIVE  DISTANCE,  X,  /xm 


Figure  4.22  Comparison  of  Theory  and  Experiments 

(Q  =  0.075W,  5o  *  5.4  nm,  k  =  5.546  x  10-3,  e  =  3.5, 
Po  =  -2.670  X  10-5.  Pi  =  1.0  X  10-2). 
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FILM  THICKNESS, 


RELATIVE  DISTANCE,  X,  yum 


Figure  4.23  Comparison  of  Theory  and  Experiments 

(Q  =  O.IW,  5o  =  4.9  nm,  k  =  5.026  x  10-3,  e  =  2.25, 
(3o  =  -2..475  X  10-4,  Pi  =  l.O  X  10-2). 
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FILM  THICKNESS 


RELATIVE  DISTANCE,  X,  /xm 


Figure  4.24  Comparison  of  Theory  and  Experiments 

(Q  =  0.125W,  5o  =  4.5  nm,  k  =  4.815  x  10-3,  e  =  1.8, 
3o  =  -6.870  X  10-4,  3j  =  1  0  X  10-2). 
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In  Table  4.1  values  of  5o.  c,  -A,  m'4,  q>d,  (T*  -  Tv).  T*  are 
given  for  the  data  presented  in  Figures  4.2-4.24.  The  values  of  -A 
were  obtained  using  Equation  (4.24.)  The  values  of  AT  =  (Ts  -  Tv) 
were  obtained  using  Equation  (4.15)  and  the  values  of  T/v  taken  to 
be  approximated  by  Ts.  The  values  of  m'4  and  q“^  were  obtained 
using: 

mid  =  aAT  (4.34) 


qid  =  aAT  Ah^  (4.35) 

We  note  that  the  values  of  AT  are  extremely  small;  too  small 
to  be  experimentally  measured.  We  also  note  the  correct  trend  in  5© 
and  q‘4  with  power  input;  as  the  power  input  increases,  8© 
decreases  and  q'4  increases.  We  further  note  that,  at  the  highest 
power  studied,  the  contact  line  thickness.  5©,  is  still  substantially 
greater  than  a  monolayer.  Higher  fluxes  were  limited  by  the 
experimental  design.  Therefore,  additional  studies  at  higher  fluxes 
with  an  improved  design  are  recommended. 


The  right-hand  side  of  Eq.  (4.21)  is  also  the  dimensionless 
mass  flux  at  the  liquid-vapor  interface,  M: 


m 


1 


mid  1  +  KT| 


(4.36) 


where  m***  =  a  (  Tj  -  Ty  )  is  the  kinetic  theory  rate  with  a  negligible 
resistance  to  conduction  in  the  liquid  (  Tj  ~  ^Iv  ^  effect  of 


disjoining  pressure  or  capillarity  on  the  vapor  pressure.  Multiplying 
the  numerator  and  the  denominator  by  the  latent  heat  of 
vaporization,  can  also  be  viewed  as  the  dimensionless  interfacial 
evaporative  heat  flux. 
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Table  4.1  Values  of  6o.  £.  -A,  m,  q'<*,  (Ts-Ty),  Ts  at  Various  Power  Inputs. 


0.50  4.0  4.0  1.124  x  10  ^2  2.81  x  10*2  4.316  x  10^  2.054  x  10-3  302.7 


Table  4.1Value£  of  5o,  e.  -A,  m.  q'**,  (Ts-Tv).  Ts  at  Various  Power  Inputs.  (Continued) 
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SECTION  5 


DISCUSSION 

5.1.  Macroscopic  Results 

Using  the  macroscopic  heat  balance  model  the  evaporative 
heat  flow  rate  based  on  two  different  sets  of  measurements  agreed. 
This  important  result  confirms  the  presence  of  evaporation  and 
gives  a  measure  of  the  contact  line  heat  sink  which  is  needed  for 
the  design  of  small  heat  exchangers.  However,  the  microscopic 
results  are  more  significant  and  are  discussed  below.  Another 
interesting  macroscopic  phenomenon  was  observed  during  the  high 
heat  flux  studies  which  are  not  reported  in  detail  herein.  The 
evaporating  meniscus  started  to  oscillate  at  higher  power  inputs  (3 
watts  and  greater).  The  oscillations  increased  with  an  increase  in 
power  input,  but  the  meniscus  maintained  an  average  shape  and 
still  supplied  cooling  liquid  to  the  hot  spot.  The  oscillations 
appeared  to  improve  the  heat  sink  capability  of  the  meniscus.  We 
note  that  this  is  a  more  complex  phenomenon  than  that  analyzed 
above  because  the  fluid  flow  rate  may  not  be  governed  by 
disjoining  and  curvature  gradients  alone.  For  a  power  input  of  4 
watts  the  temperature  of  the  interline  was  53°C.  At  this  high 
temperature,  large  surface  tension  gradients  may  start  to  influence 
the  liquid  flow  rate  [45]. 

5.2.  Microscopic  Results 

5.2.1.  Comparison  of  Data  and  Model 

By  the  proper  selection  of  e,  po  and  3 1 ,  good  agreement  could 
be  obtained  between  the  model  Equation  (4.21)  and  the  data.  There 
are  two  ways  to  evaluate  this  result.  First,  the  agreement  between 


the  experimental  and  theoretical  profiles  is  good.  This  is  an 
important  achievement.  Second,  the  values  of  the  modified 
Hamaker  constant,  A,  are  not  a  function  of  the  evaporation  rate.  We 
note  that  there  is  scatter  in  the  resulting  values  of  A.  However,  we 
feel  that  this  scatter  can  be  reduced  by  further  adjustments  in  e  and 
the  ratio  Po/Pl-  We  also  feel  that  there  is  sufficient  consistency 
between  the  systems  to  claim  success.  Although  the  polar  and  non¬ 
polar  systems  are  different,  the  subsystems  (polar  versus  non¬ 
polar)  are  consistent.  Since  this  is  the  first  nonequilibrium  use  of 
the  procedures  and  experimental  design,  direct  comparisons  with 
past  work  cannot  be  made. 

With  regard  to  the  absolute  value  of  the  Hamaker  constant,  we 
find  it  below  the  theoretical  value  of  an  ideal  system.  We  feel  that 
this  is  due  to  the  complexity  of  the  system  needed  to  study  transport 
phenomena.  Therefore,  we  have  measured  an  important  property 
(A)  in  situ  at  the  start  and  used  it  to  describe  the  characteristics  of 
an  evaporating  meniscus.  As  discussed  in  the  next  subsection,  these 
characteristics  agree  with  theoretical  expectations.  The  resulting 
confirmed  model  can  then  be  used  to  predict  the  performance  of  an 
evaporating  meniscus  as  discussed  in  Appendix  B. 

5.2.2.  Calculation  of  Dimensionless  Pressure  and 
Curvature  Profiles 

In  Figures  (5.1-5.23)  the  dimensionless  pressure,  0,  and 
curvature,  (ri"),  profiles  are  presented  for  the  various  fluids  tested. 
The  numerical  procedures  described  in  Section  4  and  Appendix  C 
were  used  to  obtain  these  results.  The  interfacial  pressure 
difference  can  be  obtained  by  multiplying  the  dimensionless 
pressure  by  the  reference  disjoining  pressure,  n©,  which  is  given  in 
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Figure  5.1  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  >  O.OW, 
8o  *  15  nm,  Jt©  *  3.76  x  10'  N/m^). 


Figure  5.2  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  »  O.IW, 
5o  *  6.4  nm,  n©  *  4.69  x  10^  N/m^). 
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Figure  5.3  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  *  0.2W, 
5o  =  6.2  nm,  Kq  =  3.65  x  10^  N/m2), 


Figure  5.4  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.3W, 
8o  =  5.9  nm,  Kq  =  5.50  x  10^  N/m^). 
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Figure  5.5  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.5W, 
So  =  5.65  nm,  jto  =  6.32  x  102  N/m2). 


Figure  5.6  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  2.0W, 
So  =  4.8  nm,  jco  «  3.62  x  102  N/m2). 
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OIMENSIONLCSS  THICKNESS 

Figure  5.7  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  O.OW, 
5o  =  11.5  nm,  iio  =  1.347  x  102  N/m^). 


DIMENSIONLESS  THICKNESS 

Figure  5.8  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.25W, 
8o  *  4,8  nm,  Xq  =  8.436  x  10^  N/m^). 
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Figure  5.9  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.35W, 
So  =  4.5  nm,  =  6.89  x  10^  N/m^). 


Figure  5.10  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.5W, 
5o  =  4.0  nm,  Jt©  »  1.756  x  10^  N/m^). 
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Figure  5.11  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  O.OW, 
So  =  13.6  nm,  Jt©  *  1-79  x  IQl  N/m2). 


Figure  5.12  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  »  0.05W, 
So  =  6.5  nm,  Ko  =  7.84  x  10*  N/m^). 
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OIMENSIONLEFS  THICKNESS 

Figure  5,13  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.07 5W, 
5o  =  5.2  nm,  “  1*22  x  10^  N/m^). 


DIMENSIONLESS  THICKNESS 

Figure  5.14  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  O.IW, 
5o  »  4.7  nm,  ti©  *  150  x  10^  N/m2). 
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Figure 


Figure 


.15  Plot  of  Dimensionless  Pressure  and  Curvatur 
vs.  Dimensionless  Film  Thickness  (Q  =  0.125V 
So  -  4.4  nm,  Jto  =  1.71  x  10^  N/m^). 


.16  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.175W, 
So  =  3.6  nm,  Jto  =  2.84  x  10^  N/m^). 
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DIMENSIONLESS  CURVATURE  DIMENSIONLESS  CURVATURE 


Figure  5.17  Plot  of  Dimensionless  Pressure  and  Curvatui 
vs.  Dimensionless  Film  Thickness  (Q  «  0.54W 
8o  =  2.9  nm,  Jto  =  5.35  x  102  N/m2). 


Figure  5.18  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  *  0.77W, 
So  ■  2.7  nm,  Jio  ■  1-85  x  lO^  N/m2). 
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DIMENSIONLESS  CURVATURE  ’  "  DIMENSIONLESS  CURVATURE 


DIMENSIONLESS  THICKNESS 


Figure  S.19  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  *  O.OW, 
So  =  17.2  nm,  iCo  =  3.88  x  lOO  N/m^). 


DIMENSIONLESS  THICKNESS 


Figure  S.20  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  »  0.05W, 
So  =  5.9  nm,  jio  *  7.53  x  10>  N/m^). 
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Figure  S.21  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  0.075W, 
5o  =  5.4  nm,  Kq  *  1.28  x  102  N/m2). 


Figure  5.22  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  =  O.IW, 
5o  *  4.9  nm,  ito  *  2.42  x  102  N/m2). 
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OtMCNSiONLCSS  THICKNESS 


Figure  5.23  Plot  of  Dimensionless  Pressure  and  Curvature 
vs.  Dimensionless  Film  Thickness  (Q  *  0.125W, 
5o  =  4.5  nm,  r©  =  3-68  x  10^  N/m^). 
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the  captions.  The  results  clearly  demonstrate  that,  initially,  fluid 
flow  results  from  a  disjoining  pressure  gradient.  In  this  region  the 
curvature  builds  up  to  a  maximum  value  at  a  relatively  small 
thickness.  The  thickness  can  be  obtained  by  multiplying  the 
dimensionless  thickness  by  the  reference  thickness,  So,  which  is  also 
given  in  each  caption.  The  effect  that  this  curvature  buildup  has  on 
the  pumping  capacity  of  a  meniscus  is  addressed  in  Appendix  B. 

5.2.3.  Calculation  of  Dimensionless  Evaporative 

Heat  or  Mass  Flux,  M,  versus  Dimensionless 
Film  Thickness 

The  dimensionless  interfacial  heat  or  mass  evaporative  flux, 

« 

M,  profiles  are  presented  in  Figures  (S.24-S.33).  The  dimensional 

evaporative  mass  flux  can  be  obtained  by  multiplying  these  values 

• 

by  the  reference  ideal  mass  flux,  m‘<^,  which  is  given  in  each  caption. 
To  obtain  the  heat  fluxes,  the  reference  ideal  heat  fluxes  are  given 
in  Table  4.1.  The  right-hand  side  of  Equation  (4.21)  which  is  equal 
to  M  ■  m/  m‘<l  «  q/q><l  can  be  used  to  describe  the  characteristics  of 
the  pronie.  At  the  interline.  So,  the  dimensionless  pressure,  is 
given  by  =  -1.  Therefore,  the  flux  is  equal  to  zero.  The  variation 
of  0  with  thickness  is  given  in  Figures  (S.l-S.23).  These  figures 
show  that  the  absolute  value  of  ^  decreases  as  the  thickness 
increases.  Therefore,  the  interfacial  flux  increases  because  the 
effect  of  interfacial  forces  decreases.  This  increase  does  not 
continue  indefinitely  because  of  conductive  resistance  in  the  liquid 
film  which  is  represented  by  kti  in  Equation  (4.21).  Therefore,  we 
obtain  a  maximum  in  the  flux  profile  and  note  that  the  resistances 
are  substantial.  All  the  trends  are  as  expected.  A  recent  numerical 
study  by  Busse,  et  al.  also  predicmd  substantial  resistances  in  the 
contact  line  region  [58]. 
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°0  50  100  150  200  250 

DIMENSIONLESS  THICKNESS 


Figure  5.24 

Dimensionless  Mass  Flux,  &f  vs.  Dimensionless  Film  Thickness,  n. 
Case  I  Q  =  O.OW,  5©  =  15  nm,  m***  =  6.83  x  lO'^  kg/m2»sec 
^Case  2  Q  =  O.IW,  8o  =  6.4  nm,  m***  =  8.50  x  10'^  kg/m2«sec^ 


Figure  5.25 

Dimensionless  Mass  Flux,  vs.  Dimensionless  Film  Thickness,  n- 
Case  I  Q  =  O.IW,  8o  “  6-^  "ni,  fn*<*  *  8.50  x  10’^  kg/m2«sec 
^Case  2  Q  a  0.3W,  8©  »  5.9  nm,  rn‘4  s  9.99  x  lO**  kg/m^«sec^ 
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®0  50  100  150  200  250 

DIMENSIONLESS  THICKNESS 


Figure  5.26 

Dimensionless  Mass  Flux,  &I  vs.  Dimensionless  Film  Thickness.  ii. 
Case  1  Q  =  0.2W,  6o  =  6.2  nm,  m**^  =  6.63  x  10*^  kg/m2»sec 
Case  2  Q  =  0.5W,  8o  =  5.65  nm,  m'd  =  1.15  x  10*3  kg/m3*sec 


DIMENSIONLESS  THICKNESS 
Figure  5.27 

Dimensionless  Mass  Flux,  vs.  Dimensionless  Film  Thickness,  n. 
I^Case  1  Q  =  O.OW,  5©  =  11.5  nm,  m**!  *  1.92  x  10'3  kg/m2»sec 
y^Case  2  Q  =  0.35W,  5©  *  4.5  nm,  m'4  a  1,04  x  lO’^  kg/m3»scc 
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Figure  5.28 

Dimensionless  Mass  Flux,  Si  vs.  Dimensionless  Film  Thickness,  n. 
Case  I  Q  =  0.25W,  8o  =  4.8  nm,  m‘4  =  1.27  x  10-2  kg/m2.sec 
^Case  2  Q  =  0.5W,  8o  =  4.0  nm,  =  2.81  x  10-2  kg/m2.sec 


®0  50  100  150  200  250 

DIMENSIONLESS  THICKNESS 


Figure  5.29 

Dimensionless  Mass  Flux,  Si  vs.  Dimensionless  Film  Thickness,  n- 
Case  1  Q  =  O.OW,  5©  =  13.6  nm,  m'4  =  5.64  x  10-6  kg/m2»sec 
Case  2  Q  »  O.IW,  5o  »  4.7  nm,  m*4  «  4.73  x  10-5  kg/m2.sec 
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Figure  5.30 

Dimensionless  Mass  Flux,  ^  vs.  Dimensionless  Film  Thickness,  t). 
Case  1  Q  =  0.075W.  So  =  5.2  nm,  mid  »  3.84  x  10*5  kg/m2«sec^ 
Case  2  Q  *  0.125W,  So  =  4.4  nm,  mid  at  5.39  x  10*5  kg/m2*sec^ 


Figure  5.31 

Dimensionless  Mass  Flux,  vs.  Dimensionless  Film  Thickness,  r\. 
Case  1  Q  »  O.OW,  So  *  17.2  nm,  mid  -  6.24  x  10*^  kg/m2»sec 
Case  2  Q  *  O.IW,  So  =  4.9  nm,  m'd  *  4.18  x  10'5  kg/m2«sec  ^ 
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®0  50  100  150  200  250 

DIMENSIONLESS  THICKNESS 


Figure  5.32 

Dimensionless  Mass  Flux.  vs.  Dimensionless  Film  Thickness,  r\. 
Case  1  Q  *  0.05W,  5o  =  5.9  nm,  m'<»  =  1.30  x  10*5  kg/m2.sec 
^  Case  2  Q  =  O.IW,  So  =  4.9  nm,  m‘4  =  4.18  x  10*5  kg/m^.sec  ^ 


°0  50  100  150  200  250 

DIMENSIONLESS  THICKNESS 

Figure  5.33 

Dimensionless  Mass  Flux,  Si  vs.  Dimensionless  Film  Thickness,  r\. 
Case  1  Q  =  0.075W,  So  »  5.4  nm,  m^d  »  2.21  x  10-5  kg/m2.sec 
^Case  2  Q  *  0.125W.  So  =  4.5  nm,  m'**  =  6.64  x  10*5  kg/m2.sec^ 
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Finally,  in  Figure  S.34,  examples  of  the  integral  evaporative  heat 
sink  for  the  same  power  input  for  the  region  studied  are  presented. 
Again,  all  the  trends  are  as  expected.  However,  the  effect  due  to 
conduction  was  not  in  the  original  model.  We  find  that  the  resistance 
due  to  thermal  conduction  also  has  a  significant  effect  on  the  integral 
evaporative  heat  sink  because  these  curves  tend  to  become  parallel  to 
the  S  axis.  Therefore,  both  interfacial  forces  and  thermal  conduction 
control  the  heat  sink  characteristics  of  the  evaporating  meniscus. 

5.2.4.  Application  of  the  Results 

We  find  that  the  simplest  overview  of  the  results  can  be 
obtained  using  the  right-hand  side  of  Equation  (4.21)  which  is  the 
dimensionless  interfacial  mass  (or  heat)  flux  given  by: 
m  1 


M  = 


m 


id 


1  +  Kr\ 


(1  +4) 


(5.1) 


This  demonstrates  that  a  smaller  (i.e.,  Kr\  small)  stable  meniscus  could 
have  a  high  evaporative  heat  flux.  The  reference  flux,  m  given  by 
Equation  (4.34)  would  be  large  with  moderate  T  and  AT.  On  the  other 
hand,  if  the  system  is  extremely  small,  the  value  of  4  (which  is  negative) 
could  make  the  flux  negligible.  In  addition,  viscous  forces  represented 
by  the  left-hand  side  of  Equation  (4.21)  could  make  the  system  unstable. 


We  find  that  the  experimental  results  confirm  the  proposed  model 
presented  in  Section  4  for  the  characteristics  of  an  evaporating 
meniscus.  Therefore,  the  model  could,  within  limits,  be  used  to  evaluate 
the  heat  sink  capability  of  other  systems.  Stability  at  higher  heat  fluxes 
is  an  unknown.  Smaller  systems  should  be  more  stable  and,  therefore, 
have  a  capacity  to  operate  at  higher  heat  fluxes.  Obviously,  this  should 
be  evaluated  with  additional  experimental  studies.  The  heat  sink 
capacity  of  transient  systems  should  also  be  evaluated  using  the  optical 
techniques  developed  within  this  study. 
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SECTION  6 


CONCLUSIONS 

1.  A  reversible  capillary  fed  heat  transfer  cell  was  successfully 
used  to  obtain  both  the  macroscopic  and  microscopic  heat 
transfer  characteristics  of  an  evaporating  extended  meniscus. 

2.  Ellipsometry  and  microcomputer  enhanced  interferometry 
were  successfully  used  to  obtain  the  thickness  profile  of  the 
meniscus  as  a  function  of  evaporation  rate. 

3.  Video  image  processing  provided  a  convenient,  efficient  and 
accurate  meth^  for  evaluating  changes  in  the  thickness 
profiles  of  evaporating  thin  liquid  films. 

4.  The  experimental  results  obtained  for  the  relatively  low  heat 
fluxes  studied  generally  confirm  the  basic  model  presented. 
Therefore,  the  model  can  be  used  to  predict  the  performance 
of  other  systems  and/or  operating  conditions. 

5.  We  find  that  there  are  significant  resistances  to  heat  transfer 
in  a  small  meniscus  due  to  interfacial  forces,  viscous  stresses, 
and  thermal  conduction.  Therefore,  the  ideal  heat  flux  based 
on  kinetic  theory  cannot  be  obtained. 
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APPENDIX  A 


A.l  Interfaciai  Phenomena  Concepts 

The  following  discussion  relates  the  concept  of  disjoining 
pressure  to  more  classical  interfaciai  concepts. 

A. 1.1  Classical 

When  a  liquid  comes  into  contact  with  a  solid  substrate,  one  of 
three  conditions  must  exist  at  equilibrium:  complete  wetting,  partial 
wetting,  or  nonwetting.  Complete  wetting  occurs  when  the  liquid 
spreads  on  the  solid  substrate  to  form  an  equilibrium  Him  with  a 
uniform  thickness  and  a  contact  angle  equal  to  zero  (  6=0).  For  0°  < 
6  <  90°  the  liquid  partially  wets  the  solid. 

At  equilibrium,  the  total  interfaciai  free  energy  is  a  minimum, 
i.e.: 

OjiAjj  +  OsvAjv  +  =  minimum  (A.l) 

where  ,  o^j,  and  Ojy  are  the  interfaciai  free  energies  of  the  solid- 

vapor,  solid-liquid  and  liquid-vapor  interfaces  respectively  and 
where  Ajj  are  the  areas.  The  solid-vapor  interfaciai  free  energy  is 

used  because  the  solid  is  in  equilibrium  with  a  vapor  phase.  The 
solid  surface  free  energy  is  a  function  of  the  vapor  pressure  because 
the  vapor  can  adsorb  on  the  solid  surface.  Equation  (A.l)  can  be 
used  to  describe  the  apparent  contact  angle  formed  by  a  liquid  in 
contact  with  a  solid  substrate: 

Ojy  cos  0  =  Ogy  -  Ojj  (A. 2) 

The  modifier  apparent  is  used  to  remind  us  that  we  cannot  see  the 
real  contact  angle,  which  is  of  molecular  dimensions,  and  that  the 
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surface  tensions  were  assumed  constant  in  the  derivation  of 
Equation  (A.2).  Since  the  interfaces  are  boundaries  between  phases, 
large  stress  gradients  are  present  in  the  contact  line  region. 
Therefore,  the  resulting  interfacial  free  energies  are  a  function  of 
the  local  liquid  film  thickness.  At  the  molecular  level,  the  interfaces 
are  dynamic  with  large  anisotropic  density  gradients.  In  passing,  it 
is  noted  that  the  internal  pressure  of  a  liquid  can  be  of  the  order  of 
a  1,000  atmospheres,  whereas  it  approximately  vanishes  in  the 
vapor  at  normal  conditions.  Obviously,  the  above  equations  are 
successful  macroscopic  models  of  the  intermolecular  force  field. 
Equivalently,  we  can  relate  the  wetting  characteristics  in  terms  of 
the  final  spreading  coefficient,  S,  introduced  by  Cooper  and  Nuttal 
[54]  : 


S=Osv*o,,.o,v 


(A.3) 


Thus  the  condition  for  complete  wetting  is  S  =  0:  while  the  condition 
for  partial  wetting  and  nonwetting  is  S  <  0. 


A. 1.2 


Disjoining  Pressure 


Although  the  apparent  contact  angle,  spreading  coefficient  and 
the  surface  tensions  give  considerable  information  concerning  the 
general  wetting  characteristics  of  a  system,  the  disjoining  pressure 
concept,  n,  outlined  next  is  more  useful.  It  allows  fluid  flow 
concepts  in  an  ultra-thin  film  to  be  introduced,  and  evaJuaied 
because  it  can  be  viewed  as  an  effective  pressure  resulting  from  a 
body  force  acting  between  the  substrate  and  the  mobile  liquid  film 
[1,2].  A  good  review  of  the  Hamaker  constant  and  disjoining 
pressure  is  given  in  [SO].  Briefly,  the  disjoining  pressure  is  (minus) 
the  potential  energy  per  unit  volume  due  to  intermolecular  forces,  F, 
and  is  a  function  of  the  film  thickness,  5; 
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n  (5)  =  -  F  (3) 


(A.4) 


Disjoining  describes  the  physical  process  whereby  a  completely 
spreading  liquid  naturally  tends  to  disjoin  a  solid  from  the  vapor  by 
spreading.  For  completeness,  an  equivalent  change  in  energy  per 
unit  surface  area,  E^,  with  liquid  film  thickness  could  be  given  as: 

n  =  -  dEj/dS  (A.5) 

Therefore,  for  a  completely  wetting  system,  the  "surface  energy" 
decreases  with  an  increase  in  the  film  thickness.  This  leads  to  a 
positive  disjoining  pressure,  n,  and  a  negative  potential  energy,  F. 

For  a  completely  spreading  system,  we  discuss  the  case  of  a 
solid  flat  plate  partially  immersed  at  an  angle  0  in  a  pool  of  wetting 
fluid  as  shown  in  Figure  2.1.  The  presence  of  interfacial  forces  in 
the  thin  wetting  film  and  capillary  forces  in  the  thicker  film 
stabilizes  it  against  the  hydrostatic  forces.  Therefore,  we  see  from 
hydrostatics  that  the  pressure,  F,  in  the  liquid  decreases  with  an 
increase  in  height.  For  a  completely  wetting  system  the  adsorbed 
thin  film  can  extend  for  a  very  long  distance.  The  interfacial  forces 
(F(3))  in  the  thin  Him  start  to  become  important  at  a  thickness  of 
approximately  10*^  m.  The  continuous  film  extends  until  it  is  a 
monolayer.  At  equilibrium,  the  chemical  potential  energy  per  unit 
volume  in  the  thin  flat  film  differs  from  that  of  a  bulk  liquid  by  an 
amount  often  referred  to  as  the  excess  potential.  Derjaguin  has 
studied  this  potential  energy  experimentally  and  theoretically,  and 
defined  it  in  units  of  a  pressure  as  the  disjoining  pressure.  The 
Dzyaloshinskii-Lifshitz-Pitaevskii  (  DLP  )  theory  enables  one  to 
calculate  the  force  (per  unit  area),  F(5),  from  the  film  thickness,  5, 
and  the  optical  and  thermophysical  properties  of  the  solid,  vapor, 
and  the  thin  liquid  film  [2].  We  note  that  although  the  DLP  theory  is 
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a  powerful  tool,  extensive  data  on  the  optical  properties  as  a 
function  of  frequency  are  needed.  An  expanded  discussion  of  this 
material  is  given  in  Reference  [36]. 

In  the  limit  of  a  thin  film  (nonretarded  regime,  5  <  lOnm)  one 
can  represent  the  dispersion  force  as; 

F(5)  =  A/  61153  ^  A  (A  6) 

while  for  the  thicker  regime  (retarded  regime,  5  >  20nm), 

F(5)  =  B/  5^  (A.7) 

In  Eq.  (A.6),  A  is  the  well  known  Hamaker  constant,  while  B  is 
defined  as  the  retarded  dispersion  force  constant.  For  stable  wetting 
films,  both  A  and  B  must  be  negative.  This  implies  that  the  potential 
energy  of  the  system  increases  when  the  wetting  film  increases  in 
thickness.  In  a  horizontal  film  or  under  zero  gravity  conditions, 
fluid  naturally  flows  from  the  thicker  to  the  thinner  region  and  a 
film  of  uniform  thickness  forms.  Calculations  of  these  forces  indicate 
that  the  above  two  equations  are  simple  approximations  but  as  such 
are  highly  useful.  The  constants  A  and  B  are  a  weak  function  of 
thickness  in  their  respective  regimes  [36]. 

The  condition  for  equilibrium  in  the  system  presented  in  Figure 
2.1  requires  that  the  following  form  of  the  chemical  potential  (per 
unit  volume)  throughout  the  liquid  be  zero; 

-F(6)  +  o  K  -pgH  =  0  (A. 8) 

where  the  three  terms  represent  the  contributions  from  the  van  der 
Waals  force,  the  capillary  force,  and  the  hydrostatic  force.  The 
curvature  of  the  film  is  represented  by  K,  the  density  by  p,  and  the 
gravitational  force  per  unit  mass  by  g.  When  we  use  only  the 


liquid-vapor  interfacial  free  energy,  we  suppress  the  subscripts. 
The  theoretical  film  profile  obtained  using  Equation  (A.8)  can  be 
compared  to  the  measured  equilibrium  profile  to  determine  the 
imponant  interactive  constant  A  or  B.  Therefore,  a  simple 
hydrostatic  experiment  can  be  done  to  determine  the  force  field  in  a 
thin  film.  The  hydrostatic  term  gives  the  pressure  as  a  function  of 
height  and,  therefore,  the  equation  gives  the  profile  as  a  function  of 
pressure.  We  note  that  both  A  and  a  are  important.  However,  the 
use  of  the  theoretically  difficult  liquid-solid  and  solid-vapor 
interfacial  free  energies  has  been  replaced  by  the  use  of  the 
Hamaker  constant.  In  the  adsorbed  thin  film  above  the  curvature 
controlled  region,  the  pressure  can  vary  by  a  large  amount  because 
H  can  be  very  large.  In  an  evaporating  system  H  and  5  decrease 
dramatically  but  the  pressure  is  still  a  function  of  the  thickness  and 
shape.  Therefore,  large  pressure  gradients  are  possible  in 
evaporating  thin  films  with  a  change  in  shape  which  usually  occurs 
over  a  very  short  distance. 

Under  zero  gravity  conditions.  Equation  (A.8)  becomes: 

-F(5)  +  oK  =  constant  (A. 9) 

Therefore,  if  an  air  or  vapor  bubble  is  present  in  a  channel  under 
zero  gravity  conditions,  its  shape  can  be  predicted.  In  a  cylindrical 
channel,  the  bubble  would  be  tubular  with  rounded  ends  and  a 
surface  film  of  liquid  (if  completely  wetting)  would  form  on  the 
walls.  The  curvature  at  the  ends  would  be  approximately  K  =  2/R 
where  R  is  the  radius  of  the  tube  if  5  «  R.  The  film  thickness  could 
then  be  obtained  from; 

A/6jt5^  or  B/5^  =  -oK 
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(A. 10) 


At  equilibrium,  the  adsorbed  thin  film  would  be  relatively  thick. 

For  example,  taking  R  »  g.6  x  I0*^m  ,  A  =  -6.9  x  10*^*  J,  o=  2  x  10'^ 
N/m  ,  we  find  5  *  2  x  10**  m. 

A.  1.2  Evaluation  of  the  Hamaker  Constant 

In  the  previous  sections,  the  central  importance  of  the 
Hamaker  constant  to  the  analysis  of  wetting,  meniscus  shape,  and 
the  heat  sink  characteristics  of  evaporating  thin  films  have  been 
discussed.  Fortunately,  the  Hamaker  "constant”  is  a  well  known 
function  of  the  film  thickness  and  the  optical  properties  of  the  liquid 
and  substrate.  However,  the  necessary  raw  data  on  the  optical 
properties  are  only  well  known  for  simple  and  ideal  systems.  For 
example,  the  properties  are  well  known  for  nonpolar  simple  fluids 
on  quartz.  In  addition  there  has  been  extensive  discussions 
concerning  water  because  of  its  importance  even  though  it  is  an 
extremely  complicated  fluid.  Some  data  for  other  polar  liquids  like 
the  alcohols  exist.  This  is  extensively  discussed  in  the  literature 
review  section  (Section  2.1). 

In  general,  the  total  disjoining  pressure,  n(5),  in  Equation  (A.4)  is 
the  sum  of  a  van  der  Waals  component  (  an  electrostatic 

component  based  on  charge  (n^|),  and  a  structural  component  (11 
The  van  der  Waals  component  is  made  up  of  Keesom  energy  (due  to 
dipole-dipole  interactions  based  on  orientation),  Debye  energy  (due 
to  dipole  -  nonpolar  interactions  based  on  polarizability),  and  the 
London  dispersion  energy  (due  to  electrodynamic  interactions 
between  the  electrons  in  all  atoms  and  molecules).  The  structural 
component  is  presumed  negligible  herein. 

The  substrate  of  primary  importance  herein  is  a  single  crystal 
of  Silicon.  A  surface  oxide  forms  on  the  Silicon  that  protects  it  from 
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the  environment,  and  the  ideal  Silicon  substrate  is  really  Silicon  with 
a  native  surface  oxide  approximately  3  nm  thick.  Theoretically,  it 
should  be  possible  to  start  with  intermolecular  forces  (as  described 
by  the  optical  properties  of  the  materials)  and  calculate  the 
Hamaker  constant,  the  ideal  heat  transfer  coefficient  and  the  heat 
sink  cf  an  evaporating  ultra-thin  film  from  first  principles. 

However,  the  exact  composition  of  the  real  experimental  interface  is 
unknown  because  it  is  a  strong  function  of  its  history  which  is  only 
approximately  known  because  of  possible  interaction  with  water 
and  other  contaminants..  Therefore,  although  all  the  theoretical 
concepts  are  available  for  overall  direction,  experimental 
characterization  of  an  "engineering"  surface  is  still  required.  The 
experimental  section  demonstrates  this  characterization. 
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APPENDIX  B:  APPLICATION  EXAMPLE 


B.l  Application  Example 

In  Figure  (B.l),  the  profiles  of  two  menisci  located  at  the  exit 
of  a  pore  are  presented.  The  pressure  jump  across  the  base  of  the 
meniscus  with  the  smaller  radius  of  curvature  is  larger  and, 
therefore,  it  would  have  the  potential  to  draw  more  liquid  to  its 
interface.  We  take  this  to  be  the  reference  profile  for  an  ideal 
evaporating  meniscus:  0  =  0  for  the  completely  wetting  case.  The 
following  material  analyzes  the  effect  that  fluid  flow  in  the  contact 
line  region  has  on  the  pressure  drop  at  the  base  of  the  meniscus. 

Briefly,  viscous  effects  in  the  immediate  vicinity  of  the  interline 
would  increase  the  effective  (apparent)  contact  angle,0g,  because  the 

meniscus  profile  is  a  function  of  the  pressure  field  that  causes  fluid 
flow  in  the  meniscus. 


The  governing  Equation  (4.21)  is  used  with  the  simplifying 
assumption  that  ^  thereby  neglecting  the  relative  effect  of 

curvature  in  the  contact  line  region. 


d^  [t\  d^J 


(1  +  Kq) 


(Bl) 


The  details  of  the  impact  of  this  assumption  are  discussed  below. 
However,  we  note  that  applying  this  restriction  in  the  contact  line 
region  does  not  effect  the  qualitative  conclusions  concerning  the 
curvature  at  the  base  of  the  meniscus. 


The  boundary  conditions  of  interest  in  this  study  are  really  a 
far-field  condition.  The  thin  film  is  a  transition  between  a  thicker 
body  and  a  nonevaporating  thin  film  (a  film  of  thickness  So  at  the 
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Ideal  Evaporating  I  Evaporating 

Meniscus  I  Meniscus 

B«sesO  r  09>os0 


B.l  Evaporating  Meniscus  Model  Profile 


# 


contact  line  at  the  exit  of  the  pore).  Therefore  the  scaled  thickness, 
n,  satisfies  the  limit: 

n  1 


as: 


Furthermore  the  slope  satisfies  the  limit: 


0 


as: 


§  — »  OO. 


The  governing  Equation  (Bl)  is  subjected  to  a  transformation 
by  defining  the  variable; 


u  = 


(B2) 


and  find  from  Equation  (Bl)  that: 

i  A  d_uA  _  1  f 

3  d^  [u  d^j  ”  1  +  Ku  1/3 


(B3) 


This  equation  is  invariant  with  respect  to  translation  in  the 
independent  variable  therefore  following  Bluman  and  Cole  [S7]  the 
transformation: 


(B4) 


is  used  with  the  chain  rule  to  eliminate  %  from  the  equation.  The 
result  is  a  first  order  differential  equation  which  may  be  integrated. 
The  proposed  solution  involves  back  transformation  to  replace  u 
with  n,  and  the  inclusion  of  the  far>field  condition: 
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^  =  [l8(I  +  K3)/n  +  <)j  + 6(^-3  -  1) 

+  9k(1  -  tt2)  +  18k2  (11-1  -  1)]'^^  (B5) 

The  solution  is  discussed  in  more  depth  in  Section  B.5.  The  integral 
heat  sink  for  the  contact  line  region,  Q,  is  directly  related  to  this 
quantity; 


Q  —  TAhm 


-^hm  ^  W 
3v/  *  u 


(B6) 


Therefore,  ~  is  a  dimensionless  contact  line  heat  sink,  Q*: 


(B7) 


B.2  Solution  for  Meniscus  Profile  and  Heat  Sink  for  Zero  k 
and  Larger  q 


Equation  (B5)  represents  a  first  order  differential  equation 
which  we  have  integrated  for  a  special  case.  The  equation  is 
simplified  to  obtain  (B9)  in  that  k  has  been  assumed  to  be  zero  (no 
insulation  effects).  It  is  funher  simplified  by  assuming  that  q  is 
large.  Specifically  we  assume: 

11*3  «l  (B8) 


Then  Equation  (B5)  becomes; 
^  =  (18  In  11  -  6)1/2 


(B9) 


Solution  of  this  is: 

q  =  ei/3  eli2/2  (BlOa) 

We  find  the  film  slope  to  be: 
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(BlOb) 


dn  1  w 
d?  =  3  7  "P 


+ 


Equation  (Bl)  was  treated  directly  by  a  numerical  method  for  zero  k 
by  Wayner,  et  al.  [9].  The  use  of  the  above  result  will  be  discussed 
more  extensively  in  Section  B.4. 


B.3  Solution  for  Meniscus  Profile  and  Heat  Sink  for  Large  k 


The  meniscus  profile  might  be  found  by  integrating  Equations 
(B4)  and  (B5).  We  propose  an  analytical  solution  valid  for  large  k. 
Analysis  of  both  Equations  (BS)  and  (B3)  for  the  case  of  large  values 
of  K  yields; 


1  ^ 
n  d^ 


k-1/2  V2 


1_  1_ 
n  *  4 


1/2 


(BID 


Integration  yields: 


/ 

T1 

V 


\ 

y 


'23+1  /  pn+if  _L 

V  I3TI-3J  18 


) 


In  this  case,  the  integral  heat  sink  is: 


Q  = 


•AAhm 

V/ 


1  ^V/2 


(B12) 


(B13) 


The  maxitnum  integral  heat  sink  is  the  heat  absorbed  by  the  entire 
meniscus,  Qt.  This  is  found  from  Equation  (B13)  with  q  "equal” 
infinity.  In  dimensional  terms: 


Qt  = 


Ahm  ky/^  /aV/^ 

V  j  tbj 


2/3 


AT 


(B14) 


B.4  Discussion 


The  expression  for  the  integral  heat  sink,  Eq.  (B6),  implies  that 
the  slope  and  curvature  of  the  Him  at  some  position  along  the 
substrate  are  functions  of  the  amount  of  heat  absorbed  by  the 
portion  of  the  meniscus  "downstream"  of  that  position,  or  in  other 
words  the  amount  of  material  flowing  past  that  position  multiplied 
by  the  heat  of  vaporization.  The  expression  for  curvature  is 
developed  by  expanding  the  left-hand  side  of  Equation  (Bl)  and 
using  the  expression  for  the  integral  heat  sink  which  provides  the 
slope  of  the  film.  That  is; 


1  d^q 


1 


q  d^ 


1 


I  +  Kq 


1  - 


(B15) 


The  relationship  between  slope,  curvature  and  the  integral  heat  sink 
may  be  illustrated  by  a  case  study. 

Consider  the  system  of  n-octane,  silicon,  and  air.  The 
temperature  is  70®C  and  the  temperature  difference  between  the 
solid  and  vapor  phases  is  S.SdO*^)^.  The  constant  -A  is  found  from 
Truong  and  Wayner  [36]  as: 

A  =  -3.18  (10-21)  J  (B16) 

It  is  imponant  to  note  that  this  is  the  value  for  an  ideal  system 
because  these  calculations  were  done  before  the  experimental 
results  were  obtained  in  this  research  project.  To  put  the  following 
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results  in  perspective  relative  to  the  experimental  portion,  the  value 
of  e  is  used.  This  demonstrates  the  importance  of  the  parameter  e. 
For  these  temperatures  and  physical  properties,  the  value  of  e  is 
0.057.  Therefore,  these  conditions  are  substantially 
different  from  those  that  were  experimentally  studied. 
Based  on  this  information  and  the  physical  properties  of  bulk  n- 
octane;  the  vapor  pressure,  Py;  the  molar  volume,  V/;  the  enthalpy  of 
vaporization,  Ahm;  and  the  kinematic  viscosity,  v; 

Pv  =  119.03  ton- 

Vy  =  1.724  (10-4)  m3/mole 

Ahm  =  339.8  kJ/kg 

v  =  4.916  (10-7)  m2/s 

We  find: 

5o  *  45.08  A 

/  =  1962  A  (B17) 

K  =  0.041 

d5  d25 

Figure  B.2  provides  a  plot  of  ^  and  versus  the  integral  heat 

sink  Q,  which  is  given  by  Equation  (B6).  Note  the  only  factor  in  Q 
which  is  varied  is  the  dimensionless  film  thickness  q.  This  plot 
shows  that  the  slope  is  steeper  and  the  curvature  larger  in  the 
portions  of  the  meniscus  which  experience  a  greater  flow.  Curiously, 
these  curves  follow  one  another  quite  closely.  The  plot  of  curvature 
is  important  because  capillarity  is  neglected  in  this  analysis.  Figure 
B.2  indicates  that  curvature  increases  with  Q.  Therefore  this  simple 
expression  for  the  integral  heat  sink  may  be  invalid  for  the  larger 
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Figure  B.2 


Film  Slope  and  Curvature  as  a  Function  of  Heat  Sink  for  a 
Particular  Case.  Table  provides  corresponding  film  thicknesses. 
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values  of  Q.  These  values  correspond  to  the  thicker  portions  of  the 
meniscus.  The  criteria  for  neglecting  capillary  pressure  are; 

d25  A 

o  (B18) 

Taking  the  surface  tension  to  be  16.87  x  10*^  N/m  and  using 
Equations  (B5)  and  (B15)  as  well  as  the  reference  quantities,  we 
estimate  the  disjoining  pressure  to  be  about  10  times  larger  than 
the  capillary  pressure  at  a  thickness  of  q  =  1.2S.  For  convenience  k 
is  set  equal  to  zero.  This  corresponds  to  Q*  =  1.04  (or  with  AT  *  4.5K, 
So  =  1  nm,  Q*  =  1.98).  A  more  in  depth  treatment  of  capillary  effects 
and  the  region  of  the  meniscus  over  which  they  are  negligible  is 
beyond  the  scope  of  this  report.  We  note  that  curvature  in  this 
example  reaches  5"  =  700,000  m*>  in  the  region  where  5  varies 
between  45 A  and  126A.  This  demonstrates  how  the  curvature 
builds  up  near  the  contact  line.  Although  the  small  scale  of  the 
analyzed  region  precludes  quantitative  comparison  with 
experimental  data  at  the  present  time,  the  results  qualitatively 
agree  with  current  experimental  observations.  The  extrapolated  use 
of  the  basic  model  evaluated  at  lower  heat  fluxes  should  be  valid. 

The  heat  flux,  q,  in  the  evaporating  meniscus  is  of  interest  as 
well.  This  may  be  estimated  with  the  ideal  evaporation  rate,  m'<i: 

q  «  m'<l  Ahin  (B19) 

This  is  illustrated  with  more  case  studies.  Case  1  corresponds  to  the 
conditions  listed  above,  with  varying  temperature  difference,  AT. 
Case  2  is  the  n-octane-silicon-air  system  at  170®C.  The  higher 
temperature  has  a  strong  effect  on  the  vapor  pressure  of  n*octane. 
At  170®C  it  is  very  near  3  atm  whereas  at  70“C  it  is  only  0.16  atm. 


The  ideal  heat  flux  is  plotted  3  .ainst  AT  in  Figure  B.3.  The  flux  is 
very  high;  however,  it  is  not  yet  apparent  how  or  if  these  fluxes 
may  be  realized.  The  effects  of  the  various  resistances  are  indicated 
in  Equation  (1.4)  and  are  detailed  in  the  main  body  of  the  report  for 
lower  heat  fluxes.  We  find  that,  because  the  interfacial  force  field 
(represented  by  lowers  the  vapor  pressure  and  conduction  in  the 
liquid  (represented  by  xq)  lowers  the  temperature  at  the  liquid- 
vapor  interface,  the  ideal  heat  flux  based  on  kinetic  theory  and  the 
solid- vapor  interfacial  temperature  difference,  AT,  is  greater  than 
the  maximum  heat  flux  in  the  evaporating  meniscus.  Therefore,  the 
results  indicate  that  high  interfacial  mass  (heat)  fluxes  can  be 
obtained  in  small  systems  (small  xq)  if  they  are  stable.  We  note  that 
stability  should  increase  in  small  wetting  systems  because  the 
relative  importance  of  the  solid  substrate  increases.  However,  this 
complex  question  has  not  been  addressed  and  the  experimental 
system  was  relatively  large. 

A  dimensionless  plot  of  the  integral  heat  sink  versus  the  film 
thickness,  constructed  from  Equation  (BS),  is  shown  in  Figure  B.4  for 
four  values  of  x.  This  plot  shows  that  the  most  active  part  of  the 
meniscus  is  the- thinnest  paru  especi^  ^r  larger  values  of  x.  Even 
though  the  heat  flux  decreases  as  q  1,  the  area  is  large  because  q' 
-*  0.  Furthermore  the  resistance  to  heat  transfer  by  the  film  is  less 
in  the  thinnest  part.  This  figure  demonstrates  the  imponance  of  x. 
Although  a  parametric  evaluation  is  required,  we  anticipate  that  x 
should  be  relatively  small  for  most  systems. 

The  analysis  may  be  related  to  the  effective  contact  angle,  6e  . 
shown  in  Figure  B.l,  due  to  evaporation  by  constructing  an 
engineering  model  of  the  evaporating  thin  film.  The  thin  film  is 
assumed  to  consist  of  two  parts.  The  thinner  pan  is  governed  by 
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Figure  B.3 


Figure  B,4 


Ideal  Heat  Flux  Versus  Solid  to  Vapor  Temperature  Change 
at  TO^C  and  170*C. 
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Integral  Heat  Sink  Versus  Film  Thickness  for  Several 
Values  of  k. 
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disjoining  pressure  and  is  analyzed  in  this  paper.  The  thicker  part  is 
a  capillary  meniscus  which  has  no  appreciable  pressure  change. 
Although  evaporation  occurs  in  the  complete  meniscus,  evaporation 
is  assumed  to  occur  only  in  the  thinner  part.  Furthermore  the 
thinner  pan  of  the  film  is  taken  to  be  so  thin  that  insulation  effects, 
KT],  are  negligible.  Therefore,  the  solution  outlined  in  Equations  (B7) 
through  (BIO)  may  be  used.  The  two  parts  of  the  film  meet  smoothly 
so  the  slope  of  the  capillary  meniscus  very  near  the  solid  is  equal  to 
the  slope  of  the  disjoining  pressure  film  at  r|.  This  slope  is  related  to 
the  integral  heat  sink  through  Equation  (BIO).  In  this  way,  the 
effective  contact  angle,  6e,  is  related  to  the  integral  heat  sink.  We 
expect  a  transition  region  between  the  capillary  meniscus  and  the 
thinner  part  of  the  film  in  which  both  effects  are  important.  We 
expect  evaporation  in  this  transition  region  and  we  expect  the  slope 
of  the  film  to  increase  and  so  the  relationship  between  the  slope  and 
the  integral  heat  sink  is  approximate.  Although  we  do  not  know  the 
accuracy  of  the  approximation,  the  additional  insight  gained  is 
significant. 


In  particular,  the  effective  contact  angle,  0e,  is  related  to  the 
dimensional  slope  through  the  inverse  tangent: 


ee 


tan-l 


8o  W 


(B20) 


According  to  this  equation  the  effective  contact  angle  increases  with 
increases  in  the  integral  heat  sink.  The  contact  angle  is  relevant 
because  it  affects  the  radius  of  curvature  of  a  meniscus  in  a  circular 
pore.  For  a  given  pore  the  radius  of  curvature,  r,  is  larger  if  the 
effective  contact  angle  0e  ‘s  larger.  Specifically: 

r  =  d/(2  cos  Oc)  (B21) 
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where  d  is  the  diameter  of  the  pore.  In  a  passive  cooling  device  the 
driving  force  for  fluid  flow  is  the  Laplace  pressure  change  at  the 
meniscus.  A  larger  value  of  r  at  the  base  of  the  meniscus  due  to  a 
change  in  8e  translates  into  a  smaller  driving  force  for  flow. 

Therefore,  a  larger  value  of  the  integral  heat  sink  may  hinder  the 
flow  of  liquid  to  the  base  of  the  meniscus  and  towards  the  interline 
from  the  base.  We  note  the  large  slopes  and  curvatures  presented 
in  Figure  B.2.  These  are  due  to  viscous  effects  because  the  film 
analyzed  would  have  been  flat  at  equilibrium. 

For  large  values  of  tc  an  analytical  expression  for  the  meniscus 

shape  was  developed.  This  is  plotted  in  Figure  B.S.  Note  that  the 

t 

independent  variable  is  ~7=r  In  the  numerical  solution  of  Wayner, 

yjlK 

et  al.  the  independent  variable  is  Therefore,  the  effect  of  large 
resistance  in  the  film  (resistance  to  thermal  conduction)  is  to  stretch 
out  the  contact  line  region.  A  stronger  effect  induces  more 
stretching. 

For  the  limiting  case  of  a  large  ic  value,  the  maximum  integral 
heat  sink  was  found  for  a  thin  film.  This  is  given  in  Equation  (B14). 
The  maximum  integral  heat  sink  is  the  total  heat  absorbed  in  the 
contact  line  region.  In  the  derivation  of  Equation  (B14)  capillary 
effects  were  neglected.  However,  this  expression  may  still  be  valid 
because,  as  previously  discussed,  the  thinnest  portions  of  the  thin 
film  are  most  active  especially  when  film  resistivity  is  large.  It  is 
the  thinnest  portions  which  have  the  least  curvature.  On  the  other 
hand,  we  caution  that  large  k  values  might  not  be 
applicable  to  heat  pipes  and  other  high  flux  systems  of 
interest. 
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Figure  B.5  Film  Profile  for  Large  k  (k  greater  than  10,  approximately). 
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The  expression.  Equation  (B14),  shows  a  reasonable 
dependence  on  parameters.  Previous  work  has  indicated  that  a  large 
value  of  the  Hamaker  constant  promotes  evaporation,  presumably 
because  a  stronger  disjoining  pressure  pumps  more  liquid  to  the 
contact  line  region.  A  higher  heat  of  vaporization  dhm  is  correlated 
with  improved  performance.  Thermal  conductivity  has  a  similar 
effect.  This  is  not  surprising  since  the  resistance  of  the  film  to 
thermal  conduction  is  assumed  to  be  very  important  in  the  derivation 
of  Equation  (B14).  Viscosity  has  a  detrimental  effect.  Lubrication 
flow,  a  viscous  flow,  involves  large  pressure  decreases  which  hinder 
evaporation.  The  parameter  "a"  indicates  the  magnitude  of  the 
response  of  evaporation  to  the  degree  of  superheating.  The 
parameter  "b"  indicates  the  propensity  of  the  hindrance  of 
evaporation  due  to  the  disjoining  pressure.  Therefore,  the 
dependence  of  Qt  on  a  and  b  is  plausible.  Finally,  we  note  that  the 
elevation  of  the  substrate  temperature  improves  the  total  heat 
uptake.  However,  we  do  not  expect  the  large  k  solution  to  yield  a  high 
heat  flux.  For  Case  I,  AT  is  5.3  (lO'^j^C  and  tc  is  0.041,  k  is  equal  to 
10  only  if  AT  is  3.6  (10''0)°C.  The  choice  of  10  is  an  estimate  based 
on  Equation  (B31).  Such  a  film  would  be  quite  thick.  So  would  be 
about  1  pm  and  the  disjoining  pressure  would  be  virtually  non¬ 
existent.  In  conclusion,  parametric  studies  are  obviously  desirable 
using  the  equations  experimentally  evaluated  herein. 

B.5  DEVELOPMENT  OF  EQUATION  (B5) 

The  governing  equation  (Bl)  is: 

d  n  1  f  1  \ 

d^  [q  d^J  “  (1  +  Kq)  '  qBj 

The  presence  of  the  q3  term  suggests  the  transformation; 
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(B23) 


u  =  n- 


This  is  only  useful  in  the  special  case  of  k  equal  to  zero.  The 
transformed  equation  is: 


1 A  fl  In'!  _  1  (  L 

3  u  d^j  ~  1  +  Ku  1/3  u 


(B24) 


As  mentioned  in  the  text  this  is  invariant  with  respect  to  translation 
in  that  is  if  we  replace  %  with  4  such  that: 


4  =  4-c 


(B25) 


where  c  is  some  constant,  the  differential  equation  for  u(^)  is 
identical  to  (B24).  Therefore  following  Bluman  and  Cole  [57]  we 


introduce  the  transformation: 
du 


(B26) 


The  chain  rule  implies: 

A.  A  w  A 

d§~d§du"  du 


(B27) 


Therefore,  equation  (B24)  may  be  rewritten: 


3  du^^uj  l  +  Ku*' 


(B28) 


IW  _  _ L_ 

3udu[uJ”  1+Ku 


_ /  ni 

1/31^  *  u  I  u 


(B29) 


This  may  be  integrated  to  yield: 
1  /wy  ,  1  / 


i(?J-  (rrU 


(B30) 
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The  integral  is  simpler  if  u  is  eliminated  in  favor  of  n  with  equation 
(B23).  We  find: 


Xfwf  ^  f  dn  j  dn 

IS^u  j  q(l  +  Kq)  ri^(l  +  Kq) 


(B31) 


The  integrals  were  solved  by  the  method  of  "partial  fractions"  for 
the  expression  (B5).  The  far-field  condition  implies  that: 

—  0  (B30) 

u 


as: 

n  -»  1  (B31) 

This  fixes  the  constant  of  integration. 
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APPENDIX  C 


C.l  Computer  Program 

The  programs  used  to  generate  the  numerical  results 
corresponding  to  the  theory  discussed  in  Section  4  of  the  report  are 
listed  in  this  appendix.  The  NAG  (Numerical  Algorithms  Group) 
routine  is  not  listed,  of  course.  The  code  consists  of  a  main  program, 
shown  in  C.2;  a  subroutine  which  calls  the  NAG  routine,  shown  in  C.3: 
and  two  subroutines  called  by  the  NAG  routine,  one  of  which 
prepares  output,  while  the  other  calculates  derivatives. 

The  main  program  reads  the  experimental  data  which  can 
be  shipped  to  a  plotting  routine  through  a  common  block,  after 
minor  alterations  for  convenience.  The  alterations  involve  reversing 
the  direction  of  the  independent  variable  (position  along  the 
substrate)  to  match  the  convention  of  the  numerical  solution, 
chopping  away  a  portion  of  the  flat  end  of  the  meniscus  (except  for 
some  occasional  noise  the  flat  end  is  strictly  flat),  and  translating  the 
data  to  reduce  the  variation  between  data  and  theory.  Reversal  of 
the  direction  of  the  independent  variable  and  translation  of  the  data 
are  justified  by  the  nature  of  the  governing  equations.  The 
dependent  variable  X  does  not  appear  and  all  terms  have  an  even 
number  of  differentiations  relative  to 

The  main  program  also  calls  subroutine  THEORl  and  can  call  a 
plotting  routine  if  desired. 

Subroutine  THEORl  calculates  the  initial  conditions  as 
discussed  in  Section  4  (note  the  direction  of  the  independent 
variable  is  reversed),  and  calls  NAG  routine  D02CBF  which  is 
reported  to  employ  Adam's  method.  Also  shown  is  subroutine  FCN, 
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called  by  ihe  NAG  routine,  which  calculates  the  derivatives  of  the 

state  variables  from  the  state  variables  U(l),  ...  U(4)  used  by  the 

^  *  . 

NAG  routine.  The  state  variables  are  defined  to  be  n.  * 

da 

-  —  in  terms  of  Section  4. 
dA. 

Different  versions  of  subroutine  "output"  are  shown  in  C.4  and 
C.5.  Each  of  these  versions  is  used  by  the  NAG  routine  to  record  the 
values  of  the  state  variables  at  particular  locations.  The  different 
versions  may  calculate  various  derived  quantities,  such  as  the 
dimensionless  curvature. 


i 
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C.2 

C  D«c.  17,  1990  Schontoerg 

C 

C  Comparison  of  Thaory  and  Sunando's  Oaca 

C 

C  Varslon  I 

c 

C  Tha  theory  coaat  from  file  ADAMSl.  The  data  is  read 

C  from  Sunando's  output  file. 

REAL’ 3  EPS , KAP. XENO , TOL, A. B , XNOT , OELREF , XREF , XTRANS 

REAL’S  X(<00» .0(400} ,XR(400) ,OR(400) , PX ( 400) , PO 1 400) ,TX(400l ,TH(4::) 

EXTERNAL  THEORl 

EXTERNAL  PLOTl 

EXTERNAL  FCN, OUTPUT 

INTEGER  NDATA, I . J , IBEG . lEND , NPLOT . HPTS 
COMMON /OERVS /EPS , RAP 
COMMON/OirrEQ/XENO, TOL. A, B, XNOT, NPTS 
COMMON / ART / OELREF . XREF . TX . TM 
COMMON / PAAM/XTRANS 
C 

REAO(S.IO)  NOATA 
:0  FORMAT (14) 

00  20  I *1. NDATA 

REAO(4.30>  X(I).0(I) 

30  FORMAT(E12.6.2X.E12.6) 

20  CONTINUE 

00  35  I>1. NOATA 
J-  NOATA  -I  *  X 
XR(J)«  X{NOATA)  -  X(I) 

OR(J)  •  0(1) 

35  CONTINUE 

READ (5. 40)  IBEG.IENO 
40  FORMAT(214) 

n?LOT’  IEND  -  IBEG  «•  1 
C 

READ (5, 60)  XTRANS 
00  50  I •IBEG. IEND 
J«  I  *  I  -  IBEG 
PX(J)  ■  XR(I)  •  XTRANS 
PO(J)  •  OR  (I) 

50  CONTINUE 
60  FORMAT(E13.6) 

C  DO  70  J'lBEG.IEND 

C  PX(J)>  PX(J)  *  XTRANS 

C70  CONTINUE 
C 

REAO(S,80)  EPS.RAP.XEND.TOL 
READ (5. 90)  A. B. XNOT 
RCAD(9.100>  OELRXF.XREP 

C  So  I  dont  have  to  cecaXcuXata  XREF  whan  I  vary  EPS  do  it  in 
C  the  progcaa. 

XREF*XREF’OSgRT(EPS) 

80  FORMAT(4E13.6) 

90  F0RMAT(3E13.6) 

100  FORMAT(2EX3.6) 

REAO(S.llO)  NPTS 
no  FORMAT  (13) 

C 

WRITE(6.X29)  EPS 

125  FORMAT) ■  CAPILLARITY  PAR  EPS  •’.£15. 6/) 

WRITE (6. 126)  RAP 

126  FORMATC  RESISTIVITY  PAR  RAP  •’,£15.6/) 

WRITE(6,X30) 

130  FORMATC  '.lOX,'  ””  RESULTS 

4'  '. 'X'.ISX. 'UIX) ’.ISX. '0(2) '.ISX. '0(3) '.ISX,  0(4) ' // 

6'  '.'REAL  XM3X.  •  H’ .  17X.  *  DHDX*  .  15X.  '  P' ,  17X. 'OPDX' /) 

WRITE (6. 132)  TOL 

132  FORMATC  TOL  •'.£15. 6/) 

WRITE(6,140)  a. a. XNOT 

140  FORMATC  A  •'.£13. 6,'  B»’,E13.6.'  XNOT- ' . E13 . 6/ /) 

CALL  THEORl 

CALL  PLOTl (PX.PO. NPLOT) 

STOP 

END 
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C.3 


i 


1 


\ 


C  D«c.  18.  1990  j«ff  schonberg 

C 

C  Subcoutin*  to  solve  extended  Keo  pcoblem 
C 

C  I  h.h.h  p'  ) •  .  3  epe  (  1  »  p  » /(  1  *  k.h  ) 

c 

C  p  •  -i/(h.h.h)  -  h’ ■ 

C 

C  The  algorithm  uses  Adam's  method. 

C 

SUBROUTINE  THEORl 

REALMS  X.XEND.TOL.U<4) ,EPS.KAP,OX,H(110) ,XNOT 
REAL'S  A.B 
REAL'S  OUMl 
COMMON / DERVS / EPS , KAB 
COMMON/OirrEQ/XEND, TOL. A. 8 . XNOT , NPTS 
EXTERNAL  FCN, OUTPUT 
C 

COMMON/OUT/OX. 1 
C 
C 

C  Now  we  try  to  compute  the  initial  conditions- 
C 

OUMl -OSQRT (EPS/ <  KAP  >  1  )> 

U<1)»  1.0  *  A'DEXP(XN0T'DUM1'3.".5)  ♦  B'DEXP  tXNOT«3 .  5> 

U<2)-  A'D0M1'(3.".5)  'DEXP(XIIOT'DOM1'3.".5) 

4*  B'(3.".5)*DEXP(XNOT'3.".5> 

U<3)«  -l.O  ♦  3,0'(  1.0  -  D0M1"2  ) 'A'DEXP(XNOT'OUMl'3. "  . 5) 

U(4)«  3.0'(  1.0  -  DUM1"2  ) 'A'OEXP(XNOT'DUMl'3 . " .  5)  '  (3 . "  .  5) 'OUMl 
C 
C 

X>0. 

N«4 

IRELAB'l 

irAIL'O 

C 

DX'(  XEND  )/(  NPTS  »  1  ) 

I  '  1 
C 

CALL  00 2CBF ( X , XENO . N . U , TOL . IRELAB , FCN . OUTPUT . 

4  H.IFAIL) 

C 

RETURN 

END 

END 

C 

c 

c 

SUBROUTINE  FCNCT.U.F) 

REAL'S  T.U(4) ,F(4> .EPS.KAP 
COMNON/OERVS/EPS, KAP 
C 

F(l)  >U(2) 

F(2)»  -(  U(l)"(-3)  )  -  U(3) 

F(3)  •U(4) 

F(4)  "EPS'S. 0'<  I.  ♦  U(3»  )'(  U<l)"(-3>  )/<  1.  *  KAP'U/l)  ) 

4  -3.0'U(2) '0(4» /U(l) 

RETURN 

END 

SUBROUTINE  OUTPUT (XSOL,U) 

REAL'S  XSOL.U(4) ,OX,X 

REAL'S  OELRXr.XREF,TX(400) .TH(400) 

COMMOM/ART/DELRBF, XREF , TX, TH 
COMNON/OUT/OX, I 

MRITE<S,10)  XSOL.Ud)  ,U(2)  ,U(3)  ,U(4) 

TX(tl -XREF'XSOL 
TH(I)  •DELRlF'Ud) 

C  MRITICr,20)  TXdJ.THdl 
20  rORMATC  ','TX«  ',E13.6,'TH"  ’.EIS.S) 

XSOL-  DX'I 
I«I'l 

10  FORMATC  ’.EIS.S.EIS.S.SEIS.S/) 

OEL«0.9 

IF(Ud)  .LT.OEL)  STOP 
return 

END 
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C.4 


c 

c 

c 


30 

c 


c 

20 


10 


SUBROUTIMC  OUTPUT (XSOL.U) 

RKA1.*8  XSOL,U(4),OX,X 

REAL*!  Ktmv.rLUX.KAP.EPS 

REAL*8  OEU«r,XREr,TXI400)  ,TH(400) 

COMHOM/ ART/OElREr , XREP, TX, TH 

COMHOM/OUT/OX.I 

COMMON / OERVS / EPS , KXP 

MR1TE(6,10)  XSOt.Ud)  ,U(2)  ,U(3l  ,0(4) 

TOT  TO  riNO  CURVATURE,  AND  ftOX  AND  WRITE  THEM  OUT 
TO  A  DATA  PILE  TO  BE  USED  BY  A  PLOTTING  ROUTINE. 

KORV-  -0(3)  -  1./  U(l)«»3 

FLUX*  (  1.  *  0(31  )/(i.  .  kAP»0(1)  ) 

WRITE(8,30)  0(1) ,U(3) ,kuRV, FLUX 
rORMAIC  ',4E1S.6/) 

TX(I) >XREF«XSOL 
TH(I)  >0ELREF*U(1) 

WRITE(7,20)  TX(I),TH(I) 

FORMATC  •,E13.6,*IH«  ',E13.6) 

XSOL*  OX* I 

1-1*1 

FORMATC  • ,E13.6,E16.6,3E19.6/) 

OEL-0.9 

IF(0(1)  .LT.DBL)  STOP 

RETURN 

END 


C.5 


SUBROUTINE  OUTPUT(XSOL, 0) 

REAL'S  XS0L,U(4) ,DX,X,GAM 
REAL'S  OELREr,XREF,TX(400) ,7H(400) 

COMMON / ART / DELREF , XREF , TX , TH 
COMMON/ OUT /DX, I 
C 

C  H«  n««G  to  9«t  hhh'OpOx  . 

C 

GAM«(  U(l)"3  )*  0(4) 

WRITE (6, 10)  XSOL.Od)  ,0(2)  ,0(3)  ,0(4)  ,GAM 
TX(I) •XREF'XSOL 
TH(I)  ■OELREF'UdI 
C  WRITER, 20)  TX(I),TH(I) 

20  FORMATC  •,'TX»  ’,E13.6,'TH-  ’,E13.6) 

XSOL«  OX* I 

I»I*1 

10  FORMATC  ',E13.6,E16.6,4EI9.6/) 

OEL'0.9 

IF(Ud)  .LT.OSL)  STOP 

RETURN 

END 
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